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ABSTRACT 

We use a large sample of galaxies from the Two Micron All Sky Survey (2MASS) and the Sloan Digital Sky 
Survey (SDSS) to calculate galaxy luminosity and stellar mass functions in the local Universe. We estimate 
corrections for passband shifting and galaxy evolution, as well as present-day stellar mass-to-light (M/L) ratios, 
by fitting the optical-near-infrared galaxy data with simple models. Accounting for the 8% galaxy overdensity in 
the SDSS early data release region, the optical and near-infrared luminosity functions we construct for this sample 
agree with most recent literature optical and near-infrared determinations within the uncertainties. We argue that 
2MASS is biased against low surface brightness galaxies, and use SDSS plus our knowledge of stellar populations 
to estimate the 'true' /^-band luminosity function. This has a steeper faint end slope and a slightly higher overall 
luminosity density than the direct estimate. Furthermore, assuming a universally-applicable stellar initial mass 
function (IMF), we find good agreement between the stellar mass function we derive from the 2MASS/SDSS data 
and that derived by Cole et al. (2001; MNRAS, 326, 255). The faint end slope slope for the stellar mass function 
is steeper than —1.1, reflecting the low stellar M/L ratios characteristic of low-mass galaxies. We estimate an 
upper limit to the stellar mass density in the local Universe il*h = 2.0 ± 0.6 x 10~ 3 by assuming an IMF as rich 
in low-mass stars as allowed by observations of galaxy dynamics in the local Universe. The stellar mass density 
may be lower than this value if a different IMF with fewer low-mass stars is assumed. Finally, we examine type- 
dependence in the optical and near-infrared luminosity functions and the stellar mass function. In agreement with 
previous work, we find that the characteristic luminosity or mass of early-type galaxies is larger than for later 
types, and the faint end slope is steeper for later types than for earlier types. Accounting for typing uncertainties, 
we estimate that at least half, and perhaps as much as 3/4, of the stellar mass in the Universe is in early-type 
galaxies. 

As an aid to workers in the field, we present in an appendix the relationship between model stellar M/L ratios 
and colors in SDSS/2MASS passbands, an updated discussion of near-infrared stellar M/L ratio estimates, and the 
volume-corrected distribution of g and /T-band stellar M/L ratios as a function of stellar mass. 

Subject headings: galaxies: luminosity function, mass function - galaxies: general — galaxies: evolution — 
galaxies: stellar content 



1. INTRODUCTION 

The distribution of galaxy luminosities and stellar masses 
in the present-day Universe is of fundamental importance for 
studying the assembly of galaxies over cosmic time, both ob- 
servationally and theoretically (e.g., Lilly et al. 1995; Lin et al. 
1999; Cole et al. 2000; Brinchmann & Ellis 2000; Somerville, 
Primack, & Faber 2001; Wolf et al. 2003). In addition to pro- 
viding the zero redshift baseline for luminosity function (LF) 
evolution, the local LF constrains powerfully much of the im- 
portant physics affecting the assembly of baryons in dark mat- 
ter halos. For example, gas accretion and cooling dominates 
the bright end of the LF, whereas feedback and photoionization 
affect primarily fainter galaxies (e.g., Cole et al. 2000; Benson 
et al. 2002). Near-infrared (NIR) luminosities of galaxies are 
particularly useful as the mass-to-light (M/L) ratios in the NIR 
vary only by a factor of two or less across a wide range of star 
formation (SF) histories (Bell & de Jong 2001, see also the Ap- 
pendix), contrasting with a factor of ten change in M/L ratio at 
the blue end of the optical regime. Therefore, NIR luminosities 
provide a cleaner estimate of galaxy stellar masses, which are 
more robustly predicted by the theoretical models (e.g., Gard- 
ner et al. 1997; Cole et al. 2001; Kochanek et al. 2001). The 
goal of this paper is to use the NIR Two Micron All Sky Sur- 



vey (2MASS; Skrutskie et al. 1997) in conjunction with optical 
data and redshifts from the Sloan Digital Sky Survey (SDSS; 
York et al. 2000) to explore the distribution of galaxy luminosi- 
ties in the optical and NIR, and to use these data to estimate the 
distribution of stellar masses in the local Universe. 

There have been a number of recent studies that have esti- 
mated LFs and mass functions (MFs), based on a number of re- 
cent large surveys. Around the knee of the LF, which represents 
the dominant contribution to the overall luminosity density, the 
agreement between the LFs from different surveys is good. In 
the optical, luminosity densities agree at typically the 20% level 
or better, accounting for differences in filter bandpasses and me- 
dian redshift (e.g., Norberg et al. 2002; Liske et al. 2003; Blan- 
ton et al. 2003c). A similar conclusion is found for the NIR 
K-bwd (e.g., Gardner et al. 1997; Cole et al. 2001; Kochanek 
et al. 2001). There are some indications that the behavior of 
the difficult-to-measure fainter galaxies may depend on envi- 
ronment (Tully et al. 2002), although these galaxies do not exist 
in sufficient numbers to contribute significantly to the luminos- 
ity density of the local Universe (e.g., Zabludoff & Mulchaey 
2000; Trentham & Tully 2002). 

Three notable exceptions to this concordance of recent LF 
measurements are the Las Campanas Redshift Survey (LCRS; 
Lin et al. 1996), the early SDSS LF from Blanton et al. (2001), 
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and the K-band LF estimate of Huang et al. (2003). The LCRS 
estimates are relatively consistent with more recent estimates of 
the optical LFs (e.g., Blanton et al. 2003c), but because of two 
offsetting effects: (;') the neglect of evolution, which biases the 
luminosity density to higher values; and (if) the use of isophotal 
magnitudes, which biases the luminosity density back down to 
lower values (Blanton et al. 2001, 2003c). Blanton et al. (2001) 
find >50% more luminosity density in the local Universe than 
more recent SDSS or Two Degree Field Galaxy Redshift Sur- 
vey (2dFGRS; Colless et al. 2001) estimates. This offset is due 
mostly to the neglect of galaxy evolution and partially to the 
use of crude ^-corrections (Blanton et al. 2003c). The differ- 
ence between the K-b&nd LF of Huang et al. (2003) and other 
local estimates is less well-understood, but could stem from the 
neglect of evolution corrections, LF fitting uncertainties and/or 
large-scale structure (their LF estimate comes from an area of 
sky 50 times smaller than the area studied in this work; we dis- 
cuss this issue in more detail in §4.2). 

Furthermore, it is unclear if the optical and NIR LFs are mu- 
tually consistent. Cole et al. (2001) compared the optical z- 
band LF from Blanton et al. (2001) with their hybrid J/K-band 
LF, finding poor agreement. Wright (2001) finds over a fac- 
tor of two offset between extrapolations from the optical LFs 
of Blanton et al. (2001) and 2MASS -derived tf-band LFs (Cole 
et al. 2001; Kochanek et al. 2001). Given the above argument 
that luminosity densities in the optical and NIR are basically 
known to within 20%, it is unclear whether this discrepancy 
can be simply accounted for by the neglect of evolution cor- 
rections by Blanton et al. (2001), or whether, for example, this 
is an indication of gross global incompleteness in TT-band LFs. 
Furthermore, the landmark stellar MFs derived by Cole et al. 
(2001) have not been, as yet, tested systematically. 

In this paper, the first in a series of papers focusing on the 
optical and NIR properties of galaxies in the local Universe, we 
use the NIR 2MASS in conjunction with optical data and red- 
shifts from SDSS to explore in detail the LFs of galaxies over 
a factor of 6 in wavelength from the M-band (0.35^m) to the 
/f-band (2. 15/xm). We then, following the methodology of Bell 
& de Jong (2000, 2001), use the constraints on the optical-NIR 
spectral energy distributions (SEDs) in conjunction with state- 
of-the-art stellar population synthesis (SPS) models to investi- 
gate in detail the stellar MF of galaxies over a factor of 1000 
in stellar mass, assuming a universally-applicable stellar initial 
mass function (IMF). We used these stellar mass estimates in 
conjunction with a statistically-estimated cold gas masses (Hi 
and H2) to construct the cold baryonic MF in the local Universe 
and the efficiency of galaxy formation (Bell et al. 2003a, see 
also, e.g., Salucci & Persic 1999). In subsequent papers we 
will examine, e.g., the Zf-band size distribution of galaxies, the 
photometric properties of a ^-selected sample, the dust con- 
tents and SF histories of disk galaxies, and the Zf-band LF of 
bulges and disks separately, amongst other goals. 

This paper is arranged as follows. In §2, we discuss the 
data, focusing on the most important sources of error and in- 
completeness. In §3, we discuss our method for deriving k- 
corrections, evolution corrections and stellar M/L ratios. In §4, 
we construct and discuss optical and NIR LFs for our sam- 
ple of galaxies. We construct stellar MFs in §5 and discuss 
these further in §6. We summarize in §7. In the Appendix, 
we present the distribution of color-derived stellar M/L ratio 
estimates as a function of galaxy mass and fits to the color- 
M/L ratio correlations in the SDSS/2MASS passbands as aids 



to workers in the field. We assume flatter = 0.3, Q\ = 0.7, and 
Hq = 100/zkms" 1 Mpc" 1 . For estimating evolution corrections, 
we assume h = 0.7. Sections 2, 3, and the Appendix go into 
considerable detail regarding the uncertainties and stellar M/L 
ratios; thus, readers interested mainly in the results should read 
§2.1 and then skip directly to §4. 

2. THE DATA, DATA QUALITY, AND SELECTION EFFECTS 
2.1. Overview 

We use the SDSS Early Data Release (EDR; Stoughton et 
al. 2002) to provide a nearly complete 13 < r < 17.5 sam- 
ple of 22679 galaxies over 414 square degrees with accurate 
ugriz fluxes and magnitudes. We match these SDSS spectro- 
scopic sample galaxies with the 2MASS extended source cata- 
log (XSC; Jarrett et al. 2000) and point source catalog (PSC) 1 . 
To match the catalogs, we choose the closest galaxy within 2"as 
the best match (for reference, the random and systematic posi- 
tional uncertainties of 2MASS and SDSS are <200 and -50 
milliarcseconds, respectively; Pier et al. 2003). In this way, we 
have a reasonably complete 13<r<17.5 sample of galax- 
ies with redshifts, 12085 of which have a match in the 2MASS 
XSC (and therefore have ugrizK fluxes, half-light radii and con- 
centrations in r and TT-bands), 6629 of which have a match in 
the 2MASS PSC (and therefore have ugrizK fluxes, and half- 
light radii and concentration parameters in r-band), and 3965 
of which have no match in either 2MASS catalog (and thus 
have the optical data only). We choose to use only the 2MASS 
TT-band at the present time. 

A complete description of these catalogs is far beyond the 
scope of this paper (see e.g., Jarrett et al. 2000; Blanton et al. 
2001; Cole et al. 2001; Stoughton et al. 2002, for more de- 
tails). Here, we discuss the most important aspects for our 
purposes: the accuracy of the magnitudes, concentrations and 
surface brightnesses, and the completeness of the catalogs. 

2.2. Magnitude Accuracy 

An important focus of this paper is the discussion of the 
offset between optical and 2MASS TT-band LFs (e.g., Wright 
2001). Because the NIR luminosity density seemed a factor 
of two below expectations, one of the principal concerns was 
a large shortfall in either the magnitudes or numbers of galax- 
ies in 2MASS. In this section, we discuss the accuracy of the 
2MASS /T-band magnitudes in detail, and briefly summarize 
the expected accuracy of magnitudes in the other passbands. 

In many respects, TT-band data from 2MASS is the ideal tool 
for constraining galaxy LFs and the stellar MF. 2 TT-band galaxy 
luminosities are five to ten times less sensitive to dust and stellar 
population effects than optical luminosities, allowing an accu- 
rate census of stellar mass in the local Universe (e.g., Bell & de 
Jong 2001). Furthermore, 2MASS covers the entire sky homo- 
geneously, with 1 % systematic variations in zero point (Niko- 
laev et al. 2000). However, in the NIR the sky background 
is roughly a factor of 100 times brighter than the mean sur- 
face brightness of luminous galaxies, and the exposure time of 
2MASS is short (7.8 seconds with a 1.3-m telescope; Skrutskie 
et al. 1997). Thus, low surface brightness (LSB) features, such 
as LSB galaxies or the outer regions of normal galaxies, may 
be missed by 2MASS. 

'http://www.ipac.caltech.edu/2mass/releases/second/doc/ancillary/pscformat.html 
2 We note that, strictly speaking, 2MASS adopts a ^,-band that peaks at 

rather shorter wavelengths than the standard /f-band, but we we will call it 

/T-band for brevity in this paper. 
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FIG. 1. — Comparison of 2MASS AT-band Kron magnitudes with deep K- 
band magnitudes for 223 galaxies from Loveday (2000). The lower panel 
shows the difference (2MASS— Loveday) in A"-band magnitudes. The bold 
line represents the mean magnitude difference 0.01 ±0.04 (0.10 ±0.02) for 
K < 1 1 (K > 1 1). The scatter is ~ 0.2 mag. 
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FIG. 2. — Comparison of 2MASS AT-band Kron magnitudes with deep K- 
band magnitudes (transformed from //-band using the 2MASS measured H-K 
color) for 1017 galaxies from papers by Gavazzi and coworkers. The lower 
panel shows the difference (2MASS-Gavazzi) in Zf-band magnitudes. The 
bold line represents the mean magnitude difference 0.06 ± 0.02 (0. 17 ± 0.01) 
for K < 10 (K > 10). The scatter is ~ 0.2 mag. 



To test how much light 2MASS misses in the LSB outer 
parts of galaxies, we compare 2MASS ,ST-band magnitudes 
from the XSC with /f-band magnitudes from deeper imag- 
ing data. Following Cole et al. (2001), we show Loveday 's 
K-bmd Kron 3 magnitude from relatively deep data (10 min- 

3 Kron (1980) magnitudes are measured in apertures that are related to the 



utes on the Cerro Tololo International Observatory 1.5-m tele- 
scope) against 2MASS TT-band Kron magnitudes (Fig. 1). At 
K < 11, 2MASS Kron magnitudes seem quite accurate, with 
a systematic offset of 0.01±0.04 mag. At fainter magnitudes, 
2MASS Kron magnitudes underestimate the true magnitude by 
0.10±0.02 mag (scatter ~0.2 mag). Cole et al. (2001) found a 
larger offset between Second Incremental Data Release 2MASS 
/f-band and total magnitudes; since the Second Incremental 
Data Release there have been improvements to the reduction 
pipeline that have improved the quality of 2MASS /f-band 
Kron magnitudes. 

We check this offset by comparison with a larger sample 
of galaxies imaged in the //-band by Gavazzi et al. (1996a,b, 
2000) and Boselli et al. (2000) 4 . We adopt the 2MASS H-K 
color to estimate the total /f-band magnitude; the typical value 
is H-K ~ 0.25, almost independent of galaxy type. The aver- 
age offset brighter (fainter) than K = 10 mag is 0.06 (0. 17) mag, 
in the sense that Gavazzi's magnitudes are slightly brighter than 
the 2MASS Kron magnitudes (see Fig. 2). We do not adopt 
this correction in this paper owing to uncertainties in transform- 
ing H-band data into A7-band. We do, nevertheless, choose to 
adopt an offset of 0.1 mag for all galaxies (not just galaxies 
with K > 11) to better match Gavazzi's offset. We note that 
magnitudes corrected in this way will be within 0. 1 mag of to- 
tal, independent of whether one compares them to Gavazzi's or 
Loveday's total magnitudes. We adopt a 0.1 mag uncertainty in 
the correction to total /f-band fluxes, added in quadrature with 
the 2MASS random magnitude error. We tested whether the 
correction to total is a function of /f-band surface brightness, as 
one could imagine that the fraction of light lost may be larger 
for lower surface brightness galaxies. We found no correlation 
between the correction to total flux and TT-band surface bright- 
ness within the errors, supporting our use of a blanket 0.1 mag 
offset. 

We also choose to match to the 2MASS PSC. There are 
very few matches to the comparison samples: 8 from Loveday 
(2000) and 8 from the sample from Gavazzi and coworkers. 
We find mean offsets of -0.9 mag (0.3 mag RMS) and -0.8 
mag (0.3 mag RMS) for the two samples. We disregard two 
outliers (with no offset and a -3 . 1 mag offset) from the eight 
of the Gavazzi sample. We account for the large PSC offsets 
by subtracting 0.85 mag from the 2MASS PSC TT-band magni- 
tudes, and setting their TT-band errors to 0.5 mag. These mag- 
nitudes are clearly of very limited use. We use them primarily 
to constrain only roughly the ^-correction, evolution correction, 
and stellar M/L ratio estimates. In particular, our choice of K- 
band magnitude limit (extinction-corrected /f-band Kron mag- 
nitude of 13.57, with the offset included after galaxy selection) 
includes only 66 galaxies from the PSC, or just over 1% of our 
/f-band selected sample. 

Because of its high signal-to-noise, SDSS Petrosian 5 ugriz 
magnitudes are expected to be accurate to better than 0.05 mag 
in a random and systematic sense (Strauss et al. 2002; Blan- 
ton et al. 2003c). Sloan papers typically make the distinction 
between preliminary magnitudes presented by the EDR in the 
natural Sloan 2.5-m telescope system, denoted u*g*r*i*z*, and 

galaxy radius (for 2MASS, not less than 5"). 

4 This sample was used to test the circular isophotal magnitudes used by 
Kochanek et al. (2001) for their A"-band derived LF. 

5 SDSS Petrosian magnitudes are estimated within an aperture that is twice 
the radius at which the local surface brightness is 1/5 of the mean surface bright- 
ness within that radius (Strauss et al. 2002). 
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the 'true' Sloan magnitudes ugriz. We denote the EDR Pet- 
rosian magnitudes ugriz for brevity. Petrosian magnitudes of 
well-resolved early-type galaxies (with close to Z? 1 / 4 law lumi- 
nosity profiles) underestimate the total flux by ~0.1 mag be- 
cause their surface brightness profiles fall off very slowly at 
large radii (Strauss et al. 2002; Blanton et al. 2003c). In this 
paper, we crudely correct for this effect by subtracting 0. 1 mag 
from the magnitude of any galaxy with an r-band concentra- 
tion parameter of c r > 2.6 (defined in the next section). While 
simplistic, it allows us to estimate the total fluxes for early-type 
galaxies to within 0.05 mag. We adopt a magnitude error of 
0.05 mag for all galaxies, added in quadrature to the (tiny) in- 
ternal SDSS random magnitude errors. Note that we apply all 
magnitude offsets after galaxy selection. 

2.3. Concentration Parameters and Surface Brightnesses 

In this paper, we study primarily the overall luminosities and 
stellar masses of galaxies, choosing not to focus on their struc- 
tural parameters, such as concentration parameter or surface 
brightness. Nevertheless, we do use concentration parameter as 
a crude discriminant between early and late-type galaxies and 
surface brightnesses when examining the completeness of the 
galaxy samples. 

We adopt as our primary morphological classifier the r-band 
concentration parameter, c r = rgo/rso, where rgo and rso are the 
circular aperture radii within which 90% and 50% of the Pet- 
rosian flux are contained, respectively. The concentration pa- 
rameter has been extensively used within the SDSS collabora- 
tion to separate between early and late-type galaxies in a rudi- 
mentary fashion; early-type galaxies have higher c r than later 
types. This is motivated by the work of Strateva et al. (2001) 
and Shimasaku et al. (2001), who find a scattered but reason- 
able correlation between qualitative morphological classifica- 
tions and c r . Strateva et al. (2001) suggest a c r > 2.6 selection 
for early-type galaxies and this cut has been adopted by Kauff- 
mann et al. (2003a). We also adopt this criterion, primarily be- 
cause it is easily reproducible, facilitating easy comparison with 
our results by other workers. Blanton et al. (2003b) note the 
sensitivity of the the concentration parameter to seeing; more 
heavily smoothed early-type galaxies appear less concentrated 
than they would be either if they were observed with better see- 
ing or were closer. Because of this, the early type definition is 
conservative; intrinsically smaller or more distant early types 
maybe misclassified as later types owing to seeing effects. 

2MASS also gives a concentration parameter c K = r 75 /r25, 
where r 7 5 and r^ are the elliptical aperture within which 75% 
and 25% of the flux are contained. In Fig. 3, we compare the 
SDSS c r to the 2MASS c K in an effort to explore systematic 
uncertainties in the use of concentration parameters as a mor- 
phological typing tool (we compare c r to color selection later 
also in §4.4). We restrict our comparison to K < 12 galaxies, 
which have sufficient S/N to estimate ck- It is clear that there 
are systematic differences between the two definitions, which 
are manifested by zero point shifts, a non-unity slope, and a 
substantial scatter. Nevertheless, making the crude approxi- 
mation that c K ~ c r + 1 (the solid line), we can compare the 
fractions classified as early-type with both definitions. Using 
c r > 2.6, we find that 401/603 K < 12 EDR galaxies are classi- 
fied as early-type. Using c K > 3.6, we find that 439/603 galax- 
ies are classified as early-type. Furthermore, 355/401 galax- 
ies (89%±7%) of r-band classified early-types are classified as 
early type using the TT-band classification. 
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FIG. 3. — AT-band concentration parameter ck = ns/ r 25 against r-band con- 
centration parameter c r = rgo/r^. The solid line denotes ck ~ c r + 1, and the 
dashed lines denote the two rough early-type cuts at c r > 2.6 and ck > 3.6. 



We use g, r, and K-bsnA surface brightnesses only as a rough 
check on the completeness properties of the 2MASS -matched 
sample and on the /f-band LF. These surface brightnesses are 
defined to be the average surface brightness within the half- 
light radii. Since the magnitudes for SDSS and 2MASS XSC 
are accurate to at worst 20%, we expect that the half-light sur- 
face brightnesses will be accurate to < 30%, given the 20% er- 
ror in total magnitude, added in quadrature with the effect of a 
20% scale size error, which is typical of scale-size comparisons 
between different authors (see, e.g., Bell & de Jong 2000). Ac- 
curacies of this order are more than sufficient for our present 
purposes. 

2.4. Completeness 

The homogeneity and completeness of SDSS and 2MASS 
make them powerful tools for understanding the characteristics 
of galaxies in the local Universe. To construct meaningful LFs 
from these datasets, we must understand the completeness char- 
acteristics of each survey. We choose SDSS EDR spectroscopic 
sample galaxies with Galactic foreground extinction-corrected 
13 < r < 17.5 (Schlegel, Finkbeiner, & Davis 1998), following 
Stoughton et al. (2002). This galaxy sample is nearly complete, 
as discussed in much more detail by e.g. Blanton et al. (2001), 
Stoughton et al. (2002) or Strauss et al. (2002). We find that the 
area covered by the SDSS spectroscopic sample is 414 square 
degrees, 90% of the 462 square degrees covered by the EDR 
imaging data (Stoughton et al. 2002). This is quite consistent 
with the the statement by Stoughton et al. (2002) that only 93% 
of the spectroscopic tiles were attempted; we adopt the 3% dif- 
ference as our systematic error in determining the sky coverage 
of this sample. We estimate a total completeness within this 
area by querying SDSS EDR photometric catalog galaxies sat- 
isfying the spectroscopic catalog inclusion criteria as outlined 
by Strauss et al. (2002). This value is 78%, which is consistent 
with a 2% loss of galaxies due to bright stars, a > 99% redshift 
success rate, and between 80% and 90% targeting efficiency 
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(e.g., Blanton et al. 2001; Stoughton et al. 2002). A value of 
85% was recently found by Nakamura et al. (2003) for bright 
SDSS galaxies (r < 15.9) in the EDR: we adopt the difference 
between our and Nakamura et al.'s measurements as the system- 
atic error in the completeness, which is propagated through into 
the <fi* and j estimates later. We do not take account of the de- 
tailed, position-dependent completeness of the sample. While 
a detailed accounting for the completeness as a function of po- 
sition is pivotal for estimating galaxy clustering properties, it 
is of only minor importance for estimating the overall LF. Fi- 
nally, we note that there is little systematic bias within SDSS 
against galaxies within the selection limits. Of order 0.1% of 
the lowest surface brightness galaxies are not targeted because 
a spectrum would be impossible to obtain, and because over 3/4 
of the lowest surface brightness features in the SDSS imaging 
survey are artifacts (Strauss et al. 2002). Also, < 5% of bright 
galaxies are rejected because they overlap a bright, saturated 
star, or because they have a very bright fiber magnitude and are 
not targeted to avoid severe cross-talk between the fiber spectra. 
Neither of these biases will significantly affect our analysis. 

In addition to estimating the completeness of SDSS inter- 
nally, we determine whether the SDSS EDR area is overdense 
using the full coverage of 2MASS. We estimate overdensities 
by comparing the number of 2MASS extended sources with 
10 < K < 13.5 in the sky outside of the Galactic Plane (\b\ > 
30°) with the number of similar sources in the SDSS EDR re- 
gion. We use an area that is slightly less than the 414 square 
degrees that we calculate for the spectroscopic EDR coverage 
because we choose rectangular areas that are fully enclosed by 
the SDSS EDR boundaries. We show in Fig. 4 that the EDR is 
overdense over the entire magnitude range 10 < K < 13.5. We 
include the estimated density in the 2dFGRS region for compar- 
ison and give the number counts for each region in Table 1 . The 
SDSS EDR is 8% overdense (with a 1% Poisson uncertainty), 
and the 2dFGRS region (used by Cole et al. 2001) is 2% un- 
derdense, compared to the whole sky. Although this estimate is 
admittedly rough because 10 < K < 13.5 galaxies are a some- 
what different set of galaxies than those with 13 < r < 17.5, 
the overdensity estimate is accurate given that we compare to 
half of the entire sky (\b\ > 30°). Furthermore, our estimate is 
insensitive to Galactic foreground extinction. We account for 
the EDR region overdensity in our analysis by multiplying the 
effective survey area by 1.08 when constructing our LFs. 

We focus our study on the 2MASS matches to the 13 < 
r < 17.5 SDSS catalog. As stated earlier, out of the 22679 
13<r<17.5 galaxies in the EDR spectroscopic sample, we 
match 12085 galaxies in the 2MASS XSC and 6629 galaxies 
in the 2MASS PSC. In Fig. 5, we explore the properties of the 
2MASS matched and unmatched galaxies in more detail. In 
the upper panels, we show the distribution of galaxy g- r color 
(left) and r-band surface brightness \i r (right). The solid Ins- 



table 1 

10 < K < 13.5 Galaxy Number Counts 



Region 


N 


Area (deg 2 ) 


«gd(deg 2 ) 


XSC \b\ > 30 


363803 


20630 


17.63 


2dFGRS 


32568 


1887 


17.26 


Sloan EDR 


7078 


369.6 


19.15 



tograms show galaxies with 2MASS counterparts, the dashed 
histograms show those without. We give the TT-band surface 
brightness /j,k (left) and apparent magnitude K (right) in the 
lower panels. We estimate \i K and K for galaxies that have no 
2MASS data using the SDSS fi r and r-band apparent magni- 
tude in conjunction with the r-K color of the best-fit SED 
model (as described in §3). We test this procedure by using 
the optical data only to predict the K-banA magnitudes of the 
12085 galaxies with K-band XSC data. We find that this proce- 
dure is accurate to 0.4 mag RMS. We see that the galaxies that 
are unmatched in 2MASS are preferentially blue and LSB in 
the optical and NIR. There are 84 LSB galaxies ((fj, K ) ~ 19.1 
magarcsec" 2 ) with estimated K < 13.57, thus, there may be 
a small population of LSB galaxies missed by 2MASS. Faint, 
LSB galaxies are visible only in the very nearest parts of an ap- 
parent magnitude-limited survey (such as 2MASS), and there- 
fore carry a large weight 1/V m . dx . Therefore, this small bias 
(~ 1%) may translate into a larger bias when considering the 
LF or luminosity density. This bias would affect all published 
2MASS LFs (e.g., Cole et al. 2001; Kochanek et al. 2001), as 
well as our own. We show later that this bias affects the faint 
end of the LF, as one would expect given the surface brightness 
dependence of the LF (see, e.g., de Jong & Lacey 2000; Cross 
& Driver 2002). We also estimate the degree of incomplete- 
ness using the optical data in conjunction with our knowledge 
of stellar populations to push the TT-band LF and stellar MF 
down to lower galaxy masses. 

We select samples for estimating LFs in different passbands 
using passband-dependent magnitude limits, following Blanton 
et al. (2003c). Specifically, when constructing ugizK LFs, we 
select the magnitude limit in ugizK so that the V max for each 
galaxy is constrained by the ugizK limit for 98% of the sam- 
ple, and is defined by the r < 17.5 limit for the other 2% of 
the galaxies. Functionally, these limits are u =18.50, g =17.74, 
i =16.94, z =16.59, and K =13.57. 

3. METHODOLOGY: ^-CORRECTIONS, EVOLUTION CORRECTIONS, 
AND STELLAR M/L RATIOS 

3.1. The Method 

To estimate LFs and stellar MFs using the redshift and ugrizK 
data for the SDSS EDR galaxies, we must estimate fe-corrections 
and stellar M/L ratios. Furthermore, Blanton et al. (2003c) and 
Norberg et al. (2002) stress the need to include the effects of 
galaxy evolution. We estimate fe-corrections, evolution correc- 
tions, and galaxy stellar M/L ratios by comparing the ugrizK 
galaxy fluxes with state-of-the art stellar population synthesis 
(SPS) models. 

For each galaxy, we construct a grid of stellar populations 
with a range of metallicities and star formation histories (SFHs) 
at both the real galaxy redshift and at redshift zero. We use 
the PEGASE model (see Fioc & Rocca-Volmerange 1997, for 
a description of an earlier version of the model), choosing ten 
galaxy metallicities from 0.5% to 250% solar. The SFHs vary 
exponentially with time t: ip = [r _1 (l -e~ T "l T )~ l ]e~'l T , where V> 
is the star formation rate (SFR), t is the exponential e-folding 
time of the SFR, and Tq is the age of the galaxy (the time since 
SF commenced). The term in the square brackets is simply 
a normalization to keep the total mass of stars formed by the 
present day at one solar mass. We choose a grid of 29 t val- 
ues between (single burst) and oo (continuous), continuing 
through to -oo and then to -1 Gyr (strongly increasing to the 
present day). Our grid covers color space relatively uniformly. 




This grid is produced at both redshift zero and at the real 
galaxy redshift. The galaxy age is 12Gyr for the redshift zero 
model, and is younger for the non-zero redshift model assuming 
h = 0.7 (in essence, we choose a formation redshift of ~ 4). For 
example, this gives an age of 10.7 Gyr for a galaxy at z = 0.1. 
We least-squares fit the model galaxies at the real galaxy red- 
shift to the observed galaxy colors to choose the best model 
galaxy template. We then estimate the evolution correction, k- 
correction, and present-day stellar M/L ratio by comparing the 
non-zero redshift model with the evolved redshift zero model. 
Thus, in essence, we correct for evolution by assuming that the 
SFH indicated by the colors of the galaxy at the observed red- 
shift continues smoothly to the present day. 

To estimate stellar masses, we adopt the z = model galaxy 
M/L ratios in each passband, assuming solar absolute mag- 
nitudes of (6.41, 5.15, 4.67, 4.56, 4.53, 3.32) in ugrizK re- 
spectively, estimated using the Pegase SPS model. Those 
wishing to convert our luminosity densities into physical units 
or SDSS or 2MASS-calibrated absolute magnitudes per cubic 
Mpc can easily use the above solar absolute magnitudes for 
conversion without loss of accuracy. We estimate uncertain- 
ties in ^-corrections, evolutionary corrections, and stellar M/L 
ratio values via three methods: (;) omitting one passband at a 
time from the SED fit (the jackknife method; 6 fits); (ii) uni- 
formly weighting all data points in the fit for each galaxy (1 
fit); and (Hi) adding random magnitude offsets with sizes cor- 
responding to the magnitude error to all the galaxy photometry 
and re-doing the fits (5 times). We then compute the errors 
from the RMS difference between these 12 different fits to the 
SED and our original SED fit. Typical ^-correction and evolu- 
tion correction random errors derived in this way for the g-band 
selected sample are (0.06, 0.03, 0.03, 0.03, 0.03, 0.02) mag in 
ugrizK respectively. Typical random errors in stellar M/L ratio 
are (0.08, 0.06, 0.05, 0.04, 0.04, 0.05) dex, again in ugrizK. 

We adopt a 'diet' Salpeter (1955) stellar IMF (following Bell 
& de Jong 2001) that has the same colors and luminosity as a 
normal Salpeter IMF, but with only 70% of the mass due to a 
lower number of faint low-mass stars. This yields stellar M/L 



ratios 30% lower at a given color than a Salpeter IMF. Bell & 
de Jong (2001) show that this IMF is 'maximum disk', inas- 
much as IMFs richer in low-mass stars over-predict the rotation 
velocity of Ursa Major Cluster galaxies with K-bwd photome- 
try and well-resolved Hi rotation curves. This prescription thus 
gives the maximum possible stellar M/L ratio. Naturally, a dif- 
ferent choice of IMF allows lower M/L ratios. For example, the 
popular Kennicutt (1983) or Kroupa, Tout & Gilmore (1993) 
IMFs have roughly 30% lower M/L ratios than this IMF, and 
are thus 'submaximal'. We discuss this issue in more detail in 
§6.1. 

3.2. Comparison with Other Constraints 

Our ^-corrections and evolution corrections are quite robust. 
TT-band fe-corrections are insensitive to galaxy spectral type. In 
particular, we find k(z) ~ -2.1 ±0.3z, which is in good agree- 
ment with k(z) ~ -2.25z from Glazebrook et al. (1995). We test 
the optical ^-corrections by comparing with a simple power- 
law interpolation, including the effects of bandpass widening. 
Blanton et al. (2003a) find that this approximation is good to 
around 0. 1 mag in all passbands, but better in riz as the spec- 
tral shapes are simpler there. We find also that our optical 
and NIR fe-corrections are consistent with the simple power- 
law recipe to within 0.1 mag in all passbands. These offsets 
decrease to 0.05 mag in riz- This agreement is more than ad- 
equate, bearing in mind our 0.05-0.1 mag fe-correction errors. 
We quantify our evolution corrections by comparing the mean 
A:+evolution correction with the mean ^-correction for our g- 
band selected galaxy sample. The mean evolution corrections 
are ~ (2.3,1.6, 1.3,1.1, 1.0,0. 8)z in ugrizK, in the sense that 
galaxies are fainter at the present day, owing mostly to pas- 
sive evolution. This can be compared to the Q values derived 
by Blanton et al. (2003c), who estimate the evolution by fit- 
ting for it explicitly in their LF estimation. They find an evolu- 
tion of - (4.2 ±0.9, 2.0 ±0.5, 1.6 ±0.3, 1.6 ±0.4,0.8 ±0.3)z in 
ugriz- Therefore, we find satisfactory agreement between our 
color-based evolutionary corrections and direct estimates from 
the LF evolution by Blanton et al. (2003c), except perhaps in 
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FIG. 5. — The g — r colors, r and A"-band surface brightnesses, and A"-band apparent magnitudes of 22679 SDSS 13 < r < 17.5 galaxies with {solid) and without 
(dashed) 2MASS matches. A"-band surface brightnesses and apparent magnitudes for the 3965 galaxies without 2MASS data are estimated using the r-band derived 
quantities and the estimated r—K color of the best-fit SED model, as described in §3. The dotted line in the upper right panel shows the surface brightness limit of 
SDSS, and the dotted line in the lower right panel shows the K < 13.57 magnitude limit that we adopt for our LF analysis. 



the M-band, where the photometric and ^-correction uncertain- 
ties are largest, and our assumption of smoothly-varying SFHs 
could easily prove inadequate. Independently, Bernardi et al. 
(2003a) find evolution of - ( 1 .2, 0.9, 0.8, 0.6)z in griz for early- 
type galaxies using a similar (but totally independent) technique 
to Blanton et al. (2003c), again within < 0.05 mag of our cor- 
rections over the redshift interval of interest. 

Bell & de Jong (2001) demonstrate that for galaxies with rel- 
atively smooth SFHs, stellar M/L ratio and optical color should 
correlate quite tightly. We present a test of our stellar M/L ra- 
tio estimates in Fig. 6. Using an independent method that ac- 
counts for bursts of SF based on the strengths of the 4000A 
break and the H<5 line, Kauffmann et al. (2003a) construct stel- 
lar M/L ratios for over 120,000 SDSS galaxies. In Fig. 19 of 
that paper, they compare their M/L ratios in g-band with the 
g-r color, estimated at z = 0.1, and find a strong correlation. 
To compare to the Kauffmann et al. (2003a) correlation we es- 
timate a color correction (g-r) z=0A ~ 0.91 (g-r), assuming a 



power-law ^-correction. Moreover, we account for the IMF dif- 
ference; our 'diet' Salpeter IMF is 0.15 dex heavier at a given 
color, because of its larger number of low-mass stars, than the 
Kennicutt IMF that Kauffmann et al. (2003a) adopt. To within 
20% random scatter our multi-color method gives results con- 
sistent with their spectral method (comparing the points with 
the solid line in Fig. 6). This is particularly impressive given 
the very different methodologies and the different stellar popu- 
lation models used 6 . 

With the low scatter in the g-r versus M/L g ratio correla- 
tion it is possible to predict stellar mass to within 20% using g- 
and r-band data alone, compared to the maximum-likelihood 
SED fits of up to 6 optical/NIR passbands. Kauffmann et al. 
(2003a) find a scatter closer to 50%; this is likely due to the 

6 Lee, McCall & Richer (2003) show a comparison of their color-derived 
two-population B-band stellar M/L ratio estimates (Lee et al. 2003) with op- 
tical color in the Appendix to that paper. They find excellent agreement with 
color-derived stellar M/L ratios from Bell & de Jong (2001), showing again that 
different methodologies yield consistent estimates of stellar M/L ratio. 
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FIG. 6. — Our g-band stellar M/L ratio estimate (from the maximum- 
likelihood fit to the galaxy SED) against g — r color for the g-band selected 
sample of 1 1848 galaxies. The dashed line is the bi-weight least-squares fit to 
the data. The solid line is a rough fit to the relationship in Fig. 19 of Kauff- 
mann et al. (2003a), accounting for the 0.15 dex offset between a Kennicutt 
IMF and our diet Salpeter IMF. In addition, we transform the Kauffmann et al. 
(2003a) z = 0. 1 color to az = color using (g— r) z ~° 1 ~ 0.91 (g—r) (see text for 
full details). The arrows show the average (short arrow) and maximum (long 
arrow) effect of dust on this relationship, as discussed in §3.3. 



different methods adopted to derive stellar M/L ratios by our 
group and Kauffmann et al. (2003a). Kauffmann et al. (2003a) 
use 3" aperture spectra, plus an r—i color that is emission-line 
sensitive for dust estimation. Thus, they are sensitive to the 
aperture mismatch between the spectra and colors, and model 
mismatches between color and spectral features. In contrast, 
we minimize the residuals explicitly between our galaxy model 
colors and the observed colors; therefore, we explicitly min- 
imize the spread in the color-M/L ratio correlation with our 
method. Either way, it is clear that we can use SDSS color data 
alone, plus a redshift, to estimate the stellar mass of galaxies to 
between 20% and 50%, relative to the answer that one obtains 
using TT-band data or spectra. In particular, this allows us to use 
the SDSS data to 'fill in' areas of parameter space not covered 
as completely by 2MASS, such as blue LSB galaxies. 

We can independently check these M/L ratios using recent 
results from Bernardi et al. (2003b). They construct estimates 
of total M/L ratio (including the potentially non-negligible con- 
tribution of dark matter) using kinematic constraints, by multi- 
plying the half-light radius by the velocity dispersion squared, 
and then dividing by half the luminosity. Assuming the Hub- 
ble Space Telescope Key Project distance scale (Freedman et al. 
2001), they then compare these M/L ratios with g-r color, find- 
ing a strong correlation (the lower right-hand panel of their Fig. 
5). Correcting for their application of a 0.08 mag bluewards 
offset in g-r color, we estimate their log 10 (M/L r ) <~ -0.15 + 
0.93(g-r). Over the g-r range of interest (0.3 < g-r < 1), 
this is within 25% at the blue end and 5% of the red end of our 
maximum-disk tuned stellar population model expectation (see 
the Appendix for details): \og m (M/L r ) = -0.306+ 1.097(g-r). 
Furthermore, their total scatter (including contributions from 
observational error) is ~ 0. 15 dex, or 40% in terms of M/L ra- 



tio, in agreement with our earlier estimate of 20% to 50%. The 
agreement between these two totally independent methodolo- 
gies, each with their own sources of systematic and random er- 
ror, is astonishing; both predict roughly a factor of five change 
in stellar M/L ratio from the blue to the red end of the galaxy 
population, and both have the same absolute stellar M/L ratio 
scale. This agreement is another, powerful argument in favor 
of a color-based stellar M/L ratio of the type discussed in the 
Appendix, or by Bell & de Jong (2001). 

3.3. Systematic uncertainties 

The above prescription for estimating ^-corrections, evolu- 
tion corrections, and stellar M/L ratios assumes that the col- 
ors of a stellar population are driven by star formation history 
(SFH) and metallicity alone. What are the systematic uncer- 
tainties introduced by neglecting the effects of dust and more 
complex SFHs? Our ^-corrections are robust, inasmuch as we 
simply use a physically-motivated model to interpolate between 
the observations (a simple power-law interpolation suffices also 
to roughly 0.1 mag). Furthermore, our evolutionary corrections 
agree with independent estimates, and since they are a relatively 
small correction (< 0.2 mag typically), small errors in the evo- 
lutionary correction will not substantially affect our results. 

However, there may be a significant uncertainty in stellar 
M/L ratio estimates that is not accounted for by our prescrip- 
tion. Overall galaxy age (i.e. the time since SF started) can 
change the stellar M/L ratio at a given color in a systematic 
sense by a small but non-negligible amount, e.g. ±0.05 dex for 
an age difference of ±3 Gyr. Furthermore, it is not a priori 
clear what effects dust may have on the stellar M/L ratios. Bell 
& de Jong (2001) show that, to first order, the effects of dust 
cancel out to within 0.1-0.2 dex when estimating color-derived 
stellar M/L ratios. This cancellation occurs because the stellar 
populations and dust each predict roughly the same amount of 
reddening per unit fading in most passbands. However, the ran- 
dom uncertainties of this technique are only 20% in terms of 
M/L ratio (see above earlier in §3), so the second order differ- 
ence between the effects of dust and stellar populations could 
be significant. 

We explore the possible effects of dust on our results in a 
simple way, following Tully et al. (1998). Tully et al. (1998) 
estimate the luminosity-dependent dust content of disk galax- 
ies by minimizing the scatter in the color-magnitude relation 
(CMR) in BRIK passbands. They find that luminous galaxies 
suffer from a 1 .7 (0.3) mag dimming in their B-band (TT-band) 
flux when going from face-on to nearly edge-on, while faint 
galaxies show very little evidence for dust. We adopt a rough 
dimming of (1.6, 1.3, 0.3) mag from face-on to edge-on (> 80°) 
in grK passbands for massive disk galaxies (masses > 3 x 
10 10 /z" 2 M Q and c r < 2.6). We allow this dimming to decrease 
linearly with logarithmic mass to zero for stellar masses below 
3 x 10 8 /i~ 2 M©. We assume a simple slab model, with an opti- 
cal depth at arbitrary inclination of r(z') = t(80°)cos80°/ cosi, 
where r(80°) is the optical depth at edge-on derived from the 
above quoted difference between edge-on and face-on. Con- 
centrated, i.e. early-type, galaxies are assigned a factor of 3 less 
dust. We assume a random distribution of orientations. 

In Fig. 6, we show schematically the effect of the average 
(short arrow) and maximum possible (long arrow) dust con- 
tents, according to our admittedly ad hoc description. The ar- 
rows show the bluing of color and reduction of M/L ratio when 
dust is taken into account. It is clear that the effects of dust and 
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stellar population are mostly degenerate in agreement with Bell 
& de Jong (2001). Nevertheless, there is a slight systematic dif- 
ference between the two effects. For this dust prescription, we 
overestimate the average stellar M/L ratio in g-band by 0.06 dex 
when we fit dust-reddened colors with pure stellar populations. 

Another source of systematic uncertainty is from bursts of 
SF. Bell & de Jong (2001) find that large bursts of SF can cause 
an over-estimate of the true M/L ratio of the stellar population, 
if the stellar population is interpreted in terms of smoothly vary- 
ing SFHs. We attempt to constrain the magnitude of this error 
for our purposes using a simple model. We choose two solar 
metallicity stellar populations, one with a decreasing SFR to the 
present day (r = 4 Gyr), and one with constant SF (r = oo Gyr). 
We then apply random variations in SFR over timescales of 
10 8 yr, distributed in a log-normal fashion with a dispersion a 
of a dex, i.e. the SFR can easily change by more than an order 
of magnitude from its baseline rate. We then examine the offset 
from the g — r color and M/L g ratio correlation of these bursty 
models, compared to smooth SFH models. For both SF mod- 
els we find that SF bursts generate a ~ 25% scatter about the 
color-M/L ratio relation, and a ~ 10% offset to slightly lower 
M/L ratio at a given color. A full order of magnitude variation 
in SFH over 10 8 yr timescales is likely to be an upper limit for 
all but the strongest present-day star-bursting galaxies; there- 
fore, we demonstrate that the bias we impose by assuming such 
simplistic SFHs is < 10%. 

To summarize, the random uncertainties of color-based stel- 
lar M/L ratio estimation are ~ 20%. Systematic uncertain- 
ties from galaxy age, dust, and bursts of SF are ~ 0.1 dex, 
or ^25%. These systematic uncertainties will not cancel out 
with larger galaxy samples, and will dominate, along with stel- 
lar IMF, the error budget of the stellar mass density of the Uni- 
verse. 

4. LUMINOSITY FUNCTIONS 

4.1. LF estimation 

We estimate LFs using the simple and intuitive V /V max for- 
malism of, e.g., Felten (1977). This method has the disad- 
vantage that it is somewhat sensitive to galaxy density varia- 
tions. For example, if the near part of the survey is rather over- 
dense, where a magnitude-limited survey is most sensitive to 
low-luminosity galaxies, then the V /V max estimator will yield 
a somewhat larger number of low-luminosity galaxies than it 
should. Both the Step-Wise Maximum Likelihood (SWML) 
method of Efstathiou et al. (1988) and the parametric method 
of Sandage et al. (1979) are insensitive to density fluctuations 
of this type (although both methods are sensitive to density fluc- 
tuations when calculating the overall LF normalization). Nev- 
ertheless, both the SWML and parametric method make the as- 
sumption that the shape of the LF is independent of environ- 
ment, yet there is impressive evidence against this assumption, 
at least in the optical (De Propris et al. 2003; Hiitsi et al. 2003). 
In contrast, the V /V max method does not make this assumption. 
Furthermore, it does not make any a priori assumptions regard- 
ing the form of the LF, unlike Sandage et al. (1979) parametric 
method. Therefore, we use the V /V max method and note, that 
among others, Cole et al. (2001) find with a similar dataset that 
LFs derived using V /V m!a and SWML are identical within the 
errors. 

For an unbiased estimate of V max , we estimate the maximum 
distance that a galaxy of a given absolute magnitude would 
be visible, accounting for Galactic foreground extinction and 



k- and evolution-corrections, not including the early-type and 
TT-band-to-total corrections. Due to slight inaccuracies in the 
k- and evolution corrections, we find a few galaxies with dis- 
tances that are larger by a few observational sigma than ex- 
pected, which gives a few galaxies with V/V max > 1. Not 
including these galaxies does not affect any of the results in 
this paper. Our formal error estimates include Poisson, Monte- 
Carlo magnitude, evolution, ^-correction, and V /V max bootstrap 
uncertainties, plus random stellar M/L ratio errors for the MFs. 
There are also systematic sources of error: e.g., the ~ 25% sys- 
tematic uncertainty from dust and bursts of SF, and the ~ 5% 
systematic uncertainty in absolute magnitude calibration (see, 
e.g., Fukugita et al. 1996). We have to first order avoided uncer- 
tainties from galaxy clustering because we have renormalized 
the luminosity functions to account for the EDR's 8% over- 
density (§2.4). However, we neither sample all of the EDR 
region, nor can we securely extrapolate to different galaxy pop- 
ulations (e.g., the 13 < r < 17.5 population as opposed to the 
10<K< 13.5 population), so we attach a 10% percent system- 
atic uncertainty to the normalization <f>* and luminosity density 
j to account for clustering. We summarize the systematic error 
budget in Table 2. 

We calculate LFs using pseudo-Mgn'z/f-limited samples (where 
98% of the galaxies are limited in the passband of interest, and 
only 2% of the galaxies are limited by their r-band flux). We 
present our results in Table 3 and include some relevant com- 
parisons from the literature. We discuss two sets of LFs in more 
detail in this paper: a joint r and TT-band selected sample, and a 
joint g and r-band selected sample such that only 2% of galaxies 
are r-band limited, and 98% are limited by the magnitude limit 
in the other passband. We do not calculate the LF or stellar MF 
for magnitude bins with less than 5 galaxies. 

4.2. The K-band limited sample 

In Fig. 7, we plot V /V max versus ZT-band absolute magni- 
tude. In an unbiased sample, an average value of V/V max = 0.5 
is expected, as galaxies uniformly fill the volume. For our sam- 
ple, the average value is 0.520 ±0.004. Excluding the 66/6282 
(1%) of galaxies with V/V max > 1 gives (V /V max ) = 0.515. This 
indicates a slight tendency for galaxies to be in the more dis- 
tant half of the sample, perhaps reflecting uncertainty in the 
evolution correction or small amounts of large-scale structure. 
Nevertheless, any bias in the sample is weak; for example, 
Cole et al. (2001) find (V /V max ) ~ 0.52 for their sample of 
2MASS/2dFGRS galaxies, yet obtain excellent agreement be- 
tween V /V max and SWML estimates of the LF. 

We plot the LFs derived using the Zf-band limited sample in 
Figs. 8 and 9. In Fig. 8, we show the TT-band LF for samples 
using different Zf-band surface brightness fix cuts. All galaxies 
in our sample (solid line) have \i K < 20 mag arcsec" 2 . We fit the 
Schechter (1976) function to the V /V max data points: 

4>(L)*L=<t>*(^) ex p(-7^)§' (D 

where <j>* is the LF normalization, L* is the characteristic lu- 
minosity at the 'knee' of the LF where the form changes from 
exponential to power law, and a is the 'faint end slope'. In 
common with other work (e.g., de Jong & Lacey 2000; Cross 
& Driver 2002), we find that a fainter surface brightness limit 
increases <f>* somewhat and substantially affects a. The K- 
band luminosity density for our galaxy sample is 5.77±0.13 x 
10 8 /iL Q Mpc" 3 (formal error only). As discussed in Table 2, we 
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Table 2 
Systematic Error Budget 



Quantity 
(1) 


Error 

(2) 


Source 
(3) 


Ref. 

(4) 


Luminosity Function 




10% 


Uncertainty in exact sky coverage (3%), completeness (7%), Poisson error in normalization 
(1%), and differences between behavior of the \§ <K < 13.5 sample and our EDR sample 


§2.4, §4.1 


M* 
a 


5% 
10% 
0.1? 


Uncertainty in absolute calibration of ugrizK system 

K only: Extrapolation to total 

Optical: from departures from a Schechter function 


(1) 

§2.2 

§4.3 


j 


+0.1 
-0.6 

15% 


NIR: from strong departures from a Schechter function, and LSB galaxy incompleteness 
Optical: from <\>* and M* uncertainty 


§2.4, §4.2 
above 




+35% 
-15% 


NIR: from </>*, M* and a uncertainty 


above & §4.2 


Stellar Mass Function 


M* & p 


30% 


Dust, bursts of SF, galaxy age, and absolute calibration uncertainty 


above & §3.3 




+0% 
-60% 


Stellar IMF 


§6.1 



References. — (1) Fukugita et al. (1996) 



Note. — Column (1) describes the quantity, (2) the contribution to the systematic error budget, (3) describes the error in more detail, and (4) gives any relevant 
references (section number or literature citation). 



Table 3 

Galaxy Luminosity Function Fits 



Band 

(1) 


mam 
(2) 


Wgal 

(3) 


(V/Vmax) 

(4) 


(z) 
(5) 


P 
(6) 


M*-51og 10 /z 
(7) 


a 
(8) 


j 

(9) 


./literature 

(10) 


Ref. 


u 


18.50 


5347 


0.532±0.004 


0.055 


0.0238(8) 


-18.13(3) 


-0.95(3) 


1.5Lto«xl0« 


1.45 xlO 8 


a 


8 


17.74 


11848 


0.509±0.003 


0.070 


0.0172(5) 


-19.73(3) 


-1.03(3) 


1-57!^ xlO 8 


1.47 xlO 8 


a 


r 


17.50 


22679 


0.509±0.002 


0.096 


0.0137(7) 


-20.57(3) 


-1.07(3) 


1.80 + o^xl0 8 


1.69 xlO 8 


a 


i 


16.94 


17984 


0.508±0.002 


0.093 


0.0118(4) 


-21.00(3) 


-1.11(3) 


2.14 + ° °2 x 10 8 


2.19xl0 8 


a 


z 


16.59 


15958 


0.520±0.002 


0.092 


0.0119(4) 


-21.34(2) 


-1.06(2) 


2.75 + ° : ?3 x 10 8 


3.22xl0 8 


a 


K 


13.57 


6282 


0.520± 0.004 


0.078 


0.0143(7) 


-23.29(5) 


-0.77(4) 


5.8 + i f x 10 8 


5.9,7xl0 8 


b,c 



Note. — The passband (1), corresponding limiting magnitude (2), and number of galaxies (3). The mean V/Vmax is in (4), and the mean redshift in (5). Each 
passband LF is fit with a Schechter function and described by three parameters - the normalization <f>* in h 3 Mpc~ 3 mag -1 (6), the characteristic luminosity L* (7), and 
the faint end slope a (8). Our luminosity density estimate j (9), compared to an estimate from the literature (10), both in units of HLq Mpc~ 3 . The literature references 
are as follows: (a) Blanton et al. (2003c); (b) Cole et al. (2001); (c) Kochanek et al. (2001). The formal error estimates for quantities are given in parentheses. Table 2 
gives a complete summary of the systematic error sources, in addition to the formal errors calculated above. The optical luminosity densities give the formal error as 
the positive error bar, and the influence of the correction of early-type galaxy magnitudes by -0.1 mag in SDSS as the negative error bar. The AT-band error estimate 
includes a substantial uncertainty from 2MASS's bias against LSB galaxies. See Fig. 15 for a graphical representation of the luminosity density literature comparison. 
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FIG. 7. — V/Vmax versus ^-band absolute magnitude. The median (thick 
solid line), and upper and lower quartiles (shaded area), are shown as a func- 
tion of X'-band absolute magnitude. The average value for the whole sample is 
0.520 ± 0.004, which is reasonably consistent with the expected value of 0.5 
(thin solid line). 
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FIG. 9. — K-band LF split by morphological type. The solid line represents 
the total LF. The dotted and dashed lines represent the LF for late and early- 
type galaxies, separated using c r = 2.6. The solid grey line denotes the 2MASS 
LF of (Kochanek et al. 2001), while the grey thick dashed and dotted lines 
denote the early/late types from that paper, respectively. 
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FIG. 8. — K-band V/Vmax LFs using different surface brightness cuts. 
The black solid line with data points represents the total sample (fix < 
20 mag arcsec~ 2 ) LF. The dotted and dashed lines represent the LF for \xk < 
17magarcsec~ 2 and fiK < 18 mag arcsec~ 2 subsamples, which shows the LF 
steepening at the faint end as the surface brightness limit gets fainter. The 
thick grey dashed line denotes the predicted ^f-band LF, in the absence of se- 
lection bias (see the text for more details). The thin grey line shows the hybrid 
Schechter+power-law fit to the predicted K-band LF. The grey solid line with 
error bars denotes the 2MASS+2dFGRS LF of (Cole et al. 2001), and the thin 
solid line is the Schechter fit to our total LF, described in Table 3. For refer- 
ence, the Schechter function fit to the total Huang et al. (2003) A'-band LF is 
shown as a dash-dotted curve. 



estimate a ±15% systematic uncertainty from our extrapolation 
to total flux, absolute magnitude calibration, and sky coverage 



uncertainty. Thus, our raw estimate of TT-band luminosity den- 
sity is 5.8±0.9 x 1O 8 /zL Mpc" 3 , including the sources of ran- 
dom and systematic error. 

Earlier, we expressed concern regarding incompleteness in 
2MASS for LSB galaxies. It is interesting to use the full 
SDSS+2MASS dataset to estimate what the /T-band LF should 
look like, in the absence of selection bias. We use the g-band 
selected galaxy sample to construct a /f-band LF using 9307 
real /f-band magnitudes and 2541 synthesized /f-band magni- 
tudes (estimated to ~ 0.4 mag accuracy using ugriz as a con- 
straint). This is denoted in Fig. 8 as a thick grey dashed line. 
The agreement is excellent at the bright end; nevertheless, the 
faint end slope of the predicted /f-band LF is substantially 
steeper. A Schechter function is a poor fit to this LF owing 
to the 'kink' at Mk -51og 10 /! ~ — 21, thus we fit a power law 
to the faint end between — 2 1 < Mk - 5 log , h < - 1 8 , which has 
a slope of a ~ -1.33. This bias against faint, LSB galaxies 
affects all 2MASS -derived estimates of not just the faint end 
slope, but also the total /f-band luminosity density. Using this 
rough hybrid Schechter+power law, which has (<p*,M*,a) = 
(0.0149, -23.33, -0.88) brightwards ofM^-51og 10 /z = -21 and 
continues with power law slope -1.33 faintwards of this limit 
(the solid grey line in Fig. 8), we estimate that the total K- 
band luminosity density may be as high as 7.6 x 10 8 hL Q Mpc -3 . 
Thus, we see that 2MASS's bias against LSB galaxies may bias 
the faint end slope downwards, and the luminosity density esti- 
mates downwards by 25%. This conclusion is qualitatively and 
quantitatively consistent with a more direct assessment of light 
missed by 2MASS's relatively shallow exposures by Andreon 
(2002) 7 . 

7 There is another argument that suggests that 2MASS misses LSB galaxies. 
The luminosity-density weighted r—K color (AB-Vega) of the galaxy popula- 
tion is 2.75±0.05 (when either r or K-band luminosity weighted), and is domi- 
nated by luminous, ~ V galaxies. SDSS should not miss large numbers of LSB 
galaxies, therefore we can use the r-band luminosity density (which should be 
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How do our luminosity functions and luminosity density esti- 
mates of x 10 8 /zL Q Mpc" 3 compare with the literature? 
In Fig. 8 (solid grey line) we compare with the LF estimate 
of Cole et al. (2001) and find excellent agreement in both the 
shape of the LF and the overall normalization. Cole et al. (2001) 
find 5.9±0.9 x 10 8 hL Q Mpc" 3 but do not account for the bias 
against LSB galaxies inherent in the 2MASS data. We com- 
pare our LF with Kochanek et al. (2001) in Fig. 9 (shown as 
a solid grey line), finding that their LF is well within our er- 
ror bars. 8 They find a somewhat steeper faint end slope than 
we do, leading them to a slightly high luminosity density of 
^7±1 x 10 8 /zL Q Mpc" 3 . Therefore, accounting for all of the 
sources of error, it is clear that we are consistent with both Cole 
et al. (2001) and Kochanek et al. (2001). Our determination has 
the advantage, however, that we have considerably reduced the 
large-scale structure uncertainty by renormalizing our LF to the 
whole \b\ > 30° sky, and that we understand in detail 2MASS's 
bias against LSB galaxies. It will indeed be interesting to see 
if large, deep Zf-band surveys will converge towards the steeper 
faint end slope predicted by our analysis. 

In this context, comparison to the rather deeper survey of 
Huang et al. (2003) is particularly interesting. They have a sam- 
ple of ~ 1000 galaxies over an area of sky 50 times smaller than 
our area, and are unable to normalize their luminosity function 
to the whole \b\ > 30° sky. Thus, their luminosity function is 
highly susceptible to the effects of large-scale structure. We 
disagree with their luminosity function (see Fig. 8) and total lu- 
minosity density of <~ 12 x 10 8 IiLq Mpc" 3 . We attribute much 
of this mismatch to large-scale structure. Two other effects may 
also contribute. First, the knee of their LF is <~ 0.2 mag brighter 
than ours, which is likely caused by the uncorrected galaxy evo- 
lution in their sample. Given a median redshift of z ~ 0.15, 
a -0.12 ±0.05 mag offset is needed to correct for this evolu- 
tion. Furthermore, ignoring even modest evolution can cause 
faint end slope over-estimation with the maximum-likelihood 
SWML and STY methods (Blanton et al. 2003c). Second, a 
Schechter fit poorly represents their dataset; an improved fit 
would have a sharper knee and a shallower faint-end slope, 
giving a substantially lower luminosity density (see Fig. 2 of 
Huang et al. 2003). However, it is intriguing that with deeper 
data they find a steep faint end slope, roughly parallel with our 
predicted TT-band LF. We predict that further work will show 
that Huang et al. (2003) indeed found roughly the right faint 
end slope, but that they were adversely affected by large-scale 
structure and a small offset from ignoring evolution corrections. 

In Fig. 9, we show the LF split crudely by morphological 
type using the SDSS r-band concentration parameter (§2.3). 
Recall that Strateva et al. (2001) and Hogg et al. (2002) show 
that most concentrated c r > 2.6 galaxies are early-type (earlier 
than Sa), although c r could be affected by seeing (see Blanton 
et al. 2003b, and §2.3). We recover the classic result that the LF 
for early types has a flat or decreasing faint end slope and has a 
brighter L* than late types, which have a somewhat steeper LF 

quite complete) plus the luminosity density-weighted r— K color (which reflects 
the behavior of ~ V galaxies) to estimate the A"-band luminosity density. This 
estimate is jk ~ 6.5 X 10 8 HLq Mpc -3 , somewhat higher than the uncorrected 
value of 5.8 X 10 8 hLQ Mpc -3 . The virtue of this r — (r - A^-based estimate 
is that the faint end slope is determined by the well-constrained r-band LF, 
whereas in the direct approach we are forced to estimate the faint end slope 
from the A'-band data directly. 

8 They also do not account for 2MASS's bias against LSB galaxies but still 
find a relatively steep faint end slope to the LF, making it hard to estimate the 
effect of the LSB bias in their case. 
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FIG. 10. — V/Vmax against g-band absolute magnitude. The median {thick 
solid line), and upper and lower quartiles {shaded area), are shown as a func- 
tion of g-band absolute magnitude. The average value for the whole sample is 
0.509 ± 0.003, reasonably consistent with the expected value of 0.5 {thin solid 
line). 



(e.g., Bromley et al. 1998; Blanton et al. 2001). Our type-split 
LFs agree qualitatively with Kochanek et al. (2001), who find 
a larger L* for early types, although their overall LF is slightly 
offset from ours. We explore the role of morphological selec- 
tion in §4.4. 

4.3. The g-band limited sample 

One strength of our combined SDSS and 2MASS sample 
is that we can construct LFs in the optical ugriz passbands to 
accompany our NIR TT-band LF. In this section, we derive a 
g-band limited LF. The Schechter fits for other passbands are 
given in Table 3 for reference, and are discussed further in §4.5. 
Analogous to the K-bwd limited sample, we show the distribu- 
tion of galaxy V /V max with absolute g-band magnitude (Fig. 
10), the g-band LF derived using different surface brightness 
limits (Fig. 1 1), and the g-band LF split by morphological type 
(Fig. 12). 

The g-band limited sample has (V /V max ) = 0.509 ±0.003, 
which is slightly higher (by <~ 3cr) than the purely random dis- 
tribution expectation of 0.5. Nevertheless, as with the K sam- 
ple, this departure is small and should not affect our results at 
more than the few percent level. 9 In Fig. 11, our g-band LF 
(solid line) compares well with that of Blanton et al. (2003c), 
shown as the dash-dotted line. Furthermore, our g-band lumi- 
nosity density (see Table 3) is ~ 7% larger than Blanton et al.'s 
value. Blanton et al. did not include light lost from the low sur- 
face brightness wings of early-type galaxies in their luminosity 
density estimate, however. When we account for the differences 

'We note that there is substantial structure in the V/V m ax distribution at 
g > -16. This structure is likely caused by large scale structure, owing to the 
small V max characteristic of galaxies with faint absolute magnitudes in apparent 
magnitude-limited samples. Furthermore, this structure is the probable origin 
of the fluctuations in the g-band LF at g > -16 in Figs. 11 and 12. 
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FIG . 11. — g-band LFs using different surface brightness cuts. The solid line 
with data points represents the total sample < 25 mag arcsec~ 2 ) LF. The 
dotted and dashed lines represent the LF for /i g < 20 mag arcsec~ 2 and fi g < 
21 magarcsec -2 subsamples, which shows the LF steepening at the faint end 
as the surface brightness limit gets fainter. The dash-dot line shows the g z=0A 
LF of Blanton et al. (2003c) transformed to redshift zero assuming unchanging 
<f>* and a, and following their Table 10. The thin solid line is the Schechter fit 
to our total LF, described in Table 3. 
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FIG. 12. — g-band LF split by morphological type. The solid line represents 
the total LF. The dotted and dashed lines represent the LF for late and early- 
type galaxies, separated using c r = 2.6. 



in technique by either neglecting the correction in our own anal- 
ysis (resulting in a reduction of our luminosity density estimate 
by 4%) or by comparing to Blanton et al.'s estimated correc- 
tion (making Blanton et al.'s estimate 3% higher), the agree- 
ment between our estimate and Blanton et al.'s is well within 
the expected uncertainties. It is worth noting that, similar to the 
TT-band LF, the inclusion of lower surface brightness galaxies 
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FIG. 14. — g-band LF of color-selected galaxy types. The solid line gives 
the total sample LF as in Fig. 12. The black dashed and dotted lines show the 
LF of color-selected early and late-type galaxies. The corresponding grey lines 
show the r-band concentration parameter-selected samples from Fig. 12. 

in g-band increases <f)* slightly and gives a steeper a. 

In Fig. 12, we show the g-band LF split into early and late 
morphological types using c r = 2.6. We find that the LF for 
early-types has a brighter L* and a flatter faint end slope than 
the later types, in agreement with many other studies of the lo- 
cal Universe (e.g., Bromley et al. 1998; Blanton et al. 2001). 
We find a flat faint-end slope for the early-type LF, which dis- 
agrees with the almost log-normal distribution seen for some 
local Universe early-type samples (e.g., Blanton et al. 2001; 
Wolf et al. 2003). Some of this discrepancy is almost certainly 
caused by a different sample selection. For example, a cut at a 
relatively red constant color will produce a log-normal LF be- 
cause of the exclusion of faint, genuinely old early-type galax- 
ies that are too blue to satisfy the color cut owing to their low 
metallicity (e.g., Bower, Lucey, & Ellis 1992). In the next sec- 
tion, we show that the early-type LF selected from a magnitude- 
dependent color cut, which accounts for the CMR, is also rela- 
tively flat (see also Bell et al. 2003b). 

4.4. Color-selection of early and late types 

We explore the role of morphological selection further by 
using the broad-band colors of galaxies. We show the CMR of 
all g-band selected galaxies in the left-hand panel of Fig. 13. 
Here, one can clearly see the 'bimodality' of the color distribu- 
tion of galaxies (e.g. Strateva et al. 2001; Blanton et al. 2003b; 
Hogg et al. 2003; Bell et al. 2003b). The galaxies separate into 
coarse blue and red 'sequences'. The blue sequence has redder 
colors at brighter magnitudes, reflecting the older ages, higher 
metallicities, and greater dust content in brighter, more mas- 
sive late-type galaxies (e.g., Tully et al. 1998; Bell & de Jong 
2000). The red sequence is also redder at brighter magnitudes, 
reflecting a metallicity-magnitude relation for older stellar pop- 
ulations (e.g., Bower, Lucey, & Ellis 1992; Kodama & Arimoto 
1997). The red sequence of galaxies from the field environ- 
ment is known to contain predominantly early-type galaxies 
(Schweizer & Seitzer 1992; Hogg et al. 2003). Thus, choos- 
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FIG. 13. — The CMR of the g-band selected galaxies (left), and the g-band selected subsample of early-types (with c r > 2.6; right). Overplotted is the 'fit' to the 
CMR with the same slope as local clusters (solid line), and the color criterion for early-type selection (dotted line), which corresponds to 0.092 mag bluer than the 
CMR fit. 



ing galaxies along the red sequence is an excellent alternative 
method for selecting early-type galaxies. 

In Fig. 13, the g,r CMR of early-type galaxies has been 
marked using a solid line with a slope of -0.03. This g — r 
slope is transformed from the CMR slope in local galaxy clus- 
ters (Bower, Lucey, & Ellis 1992; Kodama & Arimoto 1997) 
of -0.08 in U — V color using the Pegase stellar population 
synthesis models. This is in good agreement with the slope 
of -0.02 to -0.04 derived by Bernardi et al. (2003c). We de- 
fine early-type galaxies as having colors redder than A(g-r) = 
-0.092 mag from CMR ridge line, i.e., everything above the 
dotted line in Fig. 13, which corresponds to A(U-V) = 0.25 
mag following Bell et al. (2003b). 

We compare our color-magnitude based definition with the 
c r > 2.6 subsample of early-type galaxies in the right-hand 
panel of Fig. 13. Clearly the majority of concentrated galaxies 
have colors that are indicative of old stellar populations. Fully 
84% of concentrated galaxies have colors that are redder than 
our color cut. Furthermore, the fraction of concentrated galax- 
ies satisfying the color cut increases towards brighter absolute 
magnitude, meaning that the overwhelming majority of the lu- 
minosity density in concentrated galaxies will be from galax- 
ies on the CMR. Conversely, 70% of the color-selected early- 
types have c r > 2.6, although concentration has limitations as a 
morphological-classifier (§2.3). 

We show the g-band LF of early and late-types defined us- 
ing our color cut in Fig. 14. Although the color selection gives 
a larger number of early-type systems, the overall differences 
between early and late-type LFs are similar to those we find us- 
ing concentration to divide our galaxy sample morphologically 
(§4.3). Thus, the basic result that late-types have fainter L* and 
steeper faint-end slopes than early-types is robust to different 
type definitions (see also Kochanek et al. 2001). 



4.5. Comparing luminosity density estimates of the local 
Universe 

In Table 3, we see that the luminosity densities we find for 
the ugri passbands agree well with the z = densities of Blanton 
et al. (2003c). They compare their LFs and luminosity densi- 
ties in detail with various other local Universe determinations, 
finding agreement to within 10% (whether or not one corrects 
for the extrapolation of early-type galaxy magnitudes to total in 
SDSS). In the z-band, we find a ~ 15% lower luminosity den- 
sity than Blanton et al. (2003c). We agree well with Blanton et 
al.'s estimate at z = 0.1 of 2.69 ±0.05 x 10 8 /zL Q Mpc" 3 ; there- 
fore, the origin of the z-band discrepancy at z ~ is likely an 
unphysically large density evolution in their analysis. Addition- 
ally, our K-band luminosity density agrees with those of Cole 
et al. (2001) and Kochanek et al. (2001). They have also com- 
pared their LFs with literature determinations and find typically 
excellent agreement. 

Wright (2001) claims an overall discrepancy between the op- 
tical LF determinations of Blanton et al. (2001) and the NIR 
LFs of Cole et al. (2001) or Kochanek et al. (2001), finding 
a factor of ~ 2 more luminosity density in the optical than in 
the NIR. Indeed, Huang et al. (2003) claim that they 'solve' 
this problem by declaring that the the local /T-band luminos- 
ity density is a factor of two higher. However, Blanton et al. 
(2003c) show that the luminosity density estimates of Blanton 
et al. (2001) were systematically too high by > 60% (more in 
the bluer passbands). We investigate this issue through a lumi- 
nosity density comparison using our galaxy sample in different 
passbands. 

We show the luminosity density of galaxies in the local Uni- 
verse in Fig. 15. The formal uncertainties are added in quadra- 
ture to a density and magnitude systematic uncertainty of 10% 
for the optical data, and 15% for the .K-band point (because 
of the additional 10% error in extrapolating to total; Table 2). 
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FIG. 15. — The luminosity density of galaxies in the local Universe. Solid 
symbols show the results from our LF determinations in the ugrizK passbands. 
Two estimates are presented in AT-band — our direct estimate (lower solid 
point) and our estimate which accounts for 2MASS's bias against LSB galaxies 
(upper solid point). We also plot a number of literature determinations of local 
luminosity density: ugriz from Blanton et al. (2003c, stars, which for u, g 
and r are under the solid circles); bj from Norberg et al. (2002, asterisk); J 
and AT-bands from Cole et al. (2001, diamonds); A"-band from Kochanek et al. 
(2001, triangle); and the A"-band estimate of Huang et al. (2003, small circle). 
We argue that the Huang et al. (2003) result is artificially high (see text). We 
overplot three stellar population models for 12 Gyr old stellar populations of 
solar metallicity, formed with exponentially-decreasing SFRs with e-folding 
times t of 2 Gyr (dotted line), 4 Gyr (solid line), and 8 Gyr (dashed line). 



The solid symbols show the results of our analysis, and the 
open symbols a variety of literature determinations of lumi- 
nosity density, including the TT-band determination of Huang 
et al. (2003), which we have argued is artificially high. Over- 
plotted are three stellar population models for 12 Gyr old stel- 
lar populations of solar metallicity, formed with exponentially- 
decreasing SFRs with e-folding times r of 2 Gyr (dotted line), 
4 Gyr (solid line) and 8 Gyr (dashed line). All of the models are 
constrained to have the same stellar mass density as we derive in 
§5. Clearly our data reproduce many of the luminosity density 
determinations in the local Universe between u and TT-bands, 
but with the dual advantages that we use one consistent dataset 
to determine the luminosity density in all passbands, and that 
we understand sources of bias in 2MASS better than previous 
work owing to the multi-passband coverage offered by SDSS 
and 2MASS. Furthermore, the shape of the cosmic mean spec- 
trum is broadly consistent with a relatively metal-rich galaxy 
with a SFH peaked at early times and decreasing to a present 
day non-zero rate. This is in excellent agreement with the work 
of Baldry et al. (2002) and Glazebrook et al. (2003). 

Fig. 15 illustrates that we have a good understanding of the 
luminosity density of the local Universe. Thus, it is worth 
spending a few words on why there has been such confusion 
surrounding its determination. Blanton et al. (2001) presented 
luminosity densities around a factor of two higher than the new, 
revised SDSS estimates; this discrepancy is caused primarily by 
evolution correction uncertainties (Blanton et al. 2003c). These 
high estimates are at variance with earlier luminosity density 
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FIG. 16. — A"-band derived stellar MF. The solid line represents the total MR 
The dotted and dashed lines represent the MF for late and early-type galaxies, 
separated using the c r = 2.6 criteria. The naked error bars denote the 2MASS 
based MF of Cole et al. (2001), corrected to our IMF. The thin solid line is our 
Schechter function fit to the MF. 



determinations at optical (see, e.g., Fig. 1 of Wright 2001) and 
NIR wavelengths (e.g., Gardner et al. 1997; Kochanek et al. 
2001). We now see that when evolution, ^-corrections, and sys- 
tematic bias are properly and carefully accounted for (as in this 
work, or Blanton et al. 2003c), the luminosity density of the lo- 
cal Universe between 3500A and 2.2^m is understood to within 
~ 20% at a given wavelength. 

5. STELLAR MASS FUNCTIONS 

One of the strengths of our combined SDSS and 2MASS 
sample is that there are accurately measured colors for all 
the sample galaxies. Under the assumption of a universally- 
applicable stellar IMF, we can then constrain the stellar M/L ra- 
tios to within 25% in a systematic sense, given the uncertainties 
in galaxy age, dust content, and the role of bursts of SF (§3.3). 
From stellar M/L ratios, we construct galaxy stellar MFs from 
both the TT-band limited and optically-limited galaxy samples, 
which allows us to explore potential sources of systematic bias 
caused by the choice of passband. 

In Fig. 16, we show an estimate of the stellar MF from 
our /f-band limited sample. The solid line shows the stellar 
MF for all galaxies. The dashed and dotted lines show early 
and late-type galaxies, respectively. A thin solid line denotes 
the Schechter function fit to the stellar MF, and the naked er- 
ror bars show the TT-band-derived stellar MF from Cole et al. 
(2001). When expressed in terms of stellar mass, the early- 
type galaxies have a higher characteristic mass M*, and a shal- 
lower faint end slope a, than the later types. Furthermore, 
the stellar MF that we derive is in excellent agreement with 
the estimate of Cole et al. (2001). Integrating under the MF, 
we find that our total TT-band stellar mass density estimate is 
5.3 ±0.1 ± 1.6 x lO 8 /zM Mpc" 3 (random and systematic er- 
rors, respectively; see §3.3), in excellent agreement with the 
-5.5 ±0.8 x 10 8 /zM Q Mpc- 3 estimate of Cole et al. (2001). 

To explore 2MASS's bias against LSB galaxies in a com- 
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Table 4 

Galaxy Stellar Mass Function Fits 



Passband 


cj>* \og l0 M*h z 


a 


P 




^Mpc-Mog^M" 1 Mq 




hM & Mpc" 3 


All Galaxies 


8 


0.0102(5) 10.70(2) 


-1.10(2) 


5.47(1 l)x 10 s 


K 


0.0133(6) 10.63(1) 


-0.86(4) 


5.26(12) xlO 8 


Early- Type Galaxies 


8 


0.0083(4) 10.62(2) 


-0.60(4) 


3.08(6)xl0 8 


8col 


0.0107(8) 10.60(4) 


-0.70(7) 


3.84(9)xl0 8 


K 


0.0089(4) 10.61(2) 


-0.52(7) 


3.19(8)xl0 8 


Late-Type Galaxies 


8 


0.0059(3) 10.51(2) 


-1.27(3) 


2.40(4) xlO 8 


8col 


0.0027(2) 10.59(3) 


-1.45(3) 


1.70(4)xl0 8 


K 


0.0071(6) 10.48(3) 


-0.94(8) 


2.10(9)xl0 8 



Note. — A 'diet' Salpeter IMF is used for the above MFs, which satisfies maximum disk constraints from galaxies in the Ursa Major galaxy cluster (Bell & de 
Jong 2001). For reference, the Stellar mass density of Cole et al. (2001), corrected to our IMF and for light missed by the 2MASS aperture, is ~5.5 X 10 8 /iMq Mpc~ 3 . 
Formal error estimates are in parentheses. These MFs are systematically uncertain at the 30% level (in terms of mass) due to uncertainties in density normalization, 
absolute magnitude calibration, SFHs, and the role of dust (§3.3). To convert these MFs to different IMFs, add a constant to the log 10 M* following §6.1. See Table 2 
for further discussion of the systematic errors. 
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Note. — These MFs are systematically uncertain at the 30% level in terms of mass due to uncertainties in the raw LFs, SFHs, and the role of dust (§3.3). To convert 
these MFs from our default 'diet' Salpeter IMF to a different stellar IMF, add a constant to the log 10 M following §6.1. 
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FIG. 17. — g-band derived stellar MF. The solid line represents the total 
MF. The black dotted and dashed lines represent the MF for late and early- 
type galaxies, separated using the c r = 2.6 criteria. The thin solid line is our 
Schechter function fit to the MF. Overplotted in grey are the K-b&nA derived 
stellar MFs for the total sample and the two morphological subsamples from 
Fig. 16. The thin black dashed and dotted lines show the g-band MFs of color- 
selected early and late-type galaxies. The data points included in this plot are 
tabulated in Table 5. 



plementary way, we show a g-band derived stellar MF in Fig. 
17. The solid, dashed, and dotted black lines show the g-band 
stellar MF for all galaxies, early-types and late-types, respec- 
tively. The grey lines show the results from Fig. 16 for the 
K-bwd limited sample. A cursory inspection shows that the 
g-band MF is relatively poorly fit by a Schechter function; the 
break in the MF is too sharp, and the stellar MF faintwards of 
Mh 2 ~3x 10 10 Mq is better fit by a single power law than a 
Schechter function. For this reason, we present the V /V max data 
points for the g-band derived stellar mass functions in Table 5. 
The g-band stellar MF shows excellent agreement at the knee of 
the MF with the /T-band stellar MF, and shows a steeper faint- 
end slope (below 10 10 /z 2 Mq), in agreement with our earlier pre- 
diction of the 'real' /f-band LF. This again argues that 2MASS 
misses faint LSB galaxies (§2.4). The g-band stellar MF con- 
tinues to much lower masses with better signal-to-noise than 
the /T-band stellar MF, showing that the stellar MF has a rela- 
tively steep faint end (a < —1.1). A steep stellar MF contrasts 
with most contemporary determinations of galaxy LFs over 
cosmologically-significant volumes, which have a > —1.1 (e.g., 
Cole et al. 2001; Norberg et al. 2002; Blanton et al. 2003c). 
We expect a steeper stellar MF because late-type galaxies have 
lower stellar M/L ratios than earlier types, and there appears to 
be an increasing contribution from later types at low luminosi- 
ties. An interesting implication of the Universal stellar MF is 
that, because of the strong variation in optical stellar M/L ra- 
tios with galaxy stellar mass (e.g., Kauffmann et al. 2003b), the 
faint end slope of optical LFs should be shallower than the faint 
end slope of NIR LFs. Deeper optical and especially NIR LFs 
will be in a good position to explore this issue in more detail. 



6. DISCUSSION 
6.1. Changing the Stellar IMF 

Bell & de Jong (2001) find that the dominant source of uncer- 
tainty in estimating stellar M/L ratios from galaxy luminosities 
and optical colors is the stellar IMF. Uncertainties in the slope 
of the stellar IMF for high-mass stars leave the color-M/L ratio 
correlation essentially unaffected. In contrast, uncertainties in 
the behavior of the stellar IMF for low-mass stars affects the 
overall M/L ratio scale. For example, an IMF richer in low- 
mass stars yields a higher stellar M/L ratio at a given color. 
Note that the uncertainties in the low-mass end of the IMF do 
not affect galaxy colors or luminosities, and hence, the color- 
M/L ratio correlation, because these stars are too faint. 

Our stellar MFs are accurate to 30%, if we understand how 
the stellar IMF behaves at low stellar mass. This accuracy 
is useful for relative comparisons such as our estimate versus 
the stellar mass density of the Universe derived by Fukugita, 
Hogan, & Peebles (1998) or Glazebrook et al. (2003). How- 
ever, to determine in an absolute sense the stellar MF, or stellar 
mass density of the Universe, or to compare with the Universal 
cold gas density, we must account for the full range of stellar 
M/L ratio uncertainty from uncertainties in stellar IMF. 

We want to quantify a reasonable range of stellar IMF un- 
certainty. Bell & de Jong (2001) place a constraint on the stel- 
lar M/L ratios by demanding that the stellar mass in the cen- 
tral parts of spiral galaxies in the Ursa Major cluster not over- 
predict their rotation velocities. This 'maximum-disk' con- 
straint forces there to be less low-mass stars than a Salpeter 
(1955) IMF, and motivates the 'diet' Salpeter IMF. Given the 
HST key project distance scale (Freedman et al. 2001), the 
Salpeter IMF is too rich in low-mass stars to satisfy dynamical 
constraints (e.g., Weiner, Sellwood, & Williams 2001; Kauff- 
mann et al. 2003a; Kranz, Slyz, & Rix 2003). This result 
is strengthened by our earlier comparison with estimates of 
the total M/L ratios of early-types by Bernardi et al. (2003b), 
who find a color-M/L ratio correlation consistent with our 
'diet' Salpeter IMF. There are considerable uncertainties in this 
kind of analysis, including aperture bias, the effects of a non- 
isotropic velocity ellipsoid on the measured velocity dispersion 
(e.g., Cretton, Rix, & de Zeeuw 2000); nevertheless, it is en- 
couraging that maximum disks and 'maximum spheroids' yield 
consistent results to the best of our knowledge. We therefore 
adopt the 'diet' Salpeter IMF as our default IMF, and the M/L 
ratios we derive from colors, assuming this IMF, will be good 
upper limits to the stellar M/L ratio. 

Of course, the stellar M/L ratio can be lower than this max- 
imal value. For example, Bottema (1993, 1997, 1999) ar- 
gue for a substantially sub-maximal M/L ratio for all disk- 
dominated galaxies based on an analysis of the vertical ve- 
locity dispersion of stars. In a similar vein, Courteau & Rix 
(1999) argue that all disks are sub-maximal, based on a lack of 
surface-brightness dependence in the luminosity-linewidth re- 
lation (Tully & Fisher 1977). In contrast, Athanassoula, Bosma, 
& Papaioannou (1987), Weiner, Sellwood, & Williams (2001), 
and Kranz, Slyz, & Rix (2003), use a variety of techniques to 
propose a scenario in which low rotation velocity galaxies are 
substantially sub-maximal (in agreement with Bottema 1997; 
Courteau & Rix 1999), and high rotation velocity (<~ mas- 
sive) galaxies are essentially maximum-disk. The latter work 
is quite consistent with a constant, maximum-disk constrained 
IMF. Thus, it is the influence of a dark matter halo that will 
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give sub-maximal M/L ratios for dark-matter dominated, LSB, 
slower rotators (e.g., de Blok & McGaugh 1998). 10 We defer a 
clearly merited more detailed discussion of this issue to a future 
paper (de Jong & Bell, in preparation). 

Attacking the problem from another angle, there is consid- 
erable and comparable uncertainty in the determinations of the 
IMF slope for low-mass stars in the Milky Way. These un- 
certainties are well-discussed in reviews by Scalo (1998) and 
Kroupa (2002). A fair assumption is that a universal IMF ex- 
ists (although see, e.g., Scalo 1998, for a differing view), and 
the range of proposed IMFs, from Salpeter (1955) to Kroupa, 
Tout & Gilmore (1993) to the 63% of maximum-disk velocity 
(Bottema 1997), should bracket the possible IMF choices. To 
convert our maximum-disk constrained 'diet' Salpeter IMF to 
Salpeter (1955), Gould, Bahcall, & Flynn (1997), Scalo (1986), 
Kroupa, Tout & Gilmore (1993), Kroupa (2002), Kennicutt 
(1983), or Bottema 63% maximal IMFs, we should add roughly 
(0.15, 0.0, -0.1, -0.15, -0.15, -0.15, -0.35) dex to the stellar 
M/L ratios predicted using the maximal IMF 1 1 . Thus, including 
the full range of systematic uncertainty, our stellar masses can 
be increased by ~ 0.1 dex, and decreased by ~ 0.45 dex, in a 
systematic sense. For comparison with other estimates of stellar 
mass density, our results are better constrained; they should be 
changed to the same IMF by adjusting the zero-point, and the 
full range of systematic uncertainty will then be ±0.1 dex. 

6.2. The Stellar Mass Density of the Universe 

It is interesting to consider at this stage the stellar mass den- 
sity of the Universe. Adopting the 'maximal' diet Salpeter IMF, 
we find that Q, t h ~ 0.00197(4) (formal error), with a ±30% sys- 
tematic error due to LF uncertainties and the effects of dust and 
bursts of SF. This estimate of il*h ~ 0.0020 ± 0.0006 is an up- 
per limit; if all galaxies are substantially sub-maximal, this es- 
timate will need to be revised downwards. Splitting these into 
early and late types (by either concentration or color-selection), 
we find that r2», Eariy /z ~ 0.0012(4) and ft* iL ate& ~ 0.0007(2). 
This estimate accounts for classification and systematic stellar 
M/L ratio uncertainties, while assuming the same stellar IMF in 
early and late types. Thus, we find that between half and three 
quarters of the stellar mass in the local Universe is in early types 
(this is robust to exactly which universal IMF is chosen, as long 
as it applies to both early and late types). Of course, much 
of the stellar mass in late-type galaxies will also be old. This 
agrees with the conclusions of Hogg et al. (2002), who find a 
very similar result, but accounting for stellar M/L ratios in a less 
elaborate way. This conclusion is in qualitative agreement with 
the shape of the 'cosmic SFH' from direct (e.g. Madau et al. 
1996; Yan et al. 1999; Blain et al. 1999), or indirect (Baldry et 
al. 2002; Glazebrooket al. 2003), estimates; all estimates agree 
that the bulk of SF in the Universe happened at early times, and 
is much slower at the present day. 

Adopting the same 'diet' Salpeter maximum-disk tuned IMF 
as we do, Cole et al. (2001) find n*h = 2.0±0.3 x 10" 3 account- 
ing for light missed by 2MASS, Persic & Salucci (1992) find 
^spirais+eiiipticais ~2x 10" 3 , Fukugita et al. (1996) find ~ 
2.6 ± 1.3 x 10" 3 , Kochanek et al. (2001) find ft*/* = 2.4 ± 
0.3 x 10" 3 , and Glazebrook et al. (2003) find fl„h values be- 

10 See Fig. 6 of Bell & de Jong (2001) for an illustration of the higher maxi- 
mum disk M/L ratio estimates for LSB galaxies, plausibly indicating dark mat- 
ter domination, even in the inner parts of the galaxies. 

"Note that the Salpeter case violates the maximum disk constraints, how- 
ever. 



tween 1.8 x 10" 3 and 3.9 x 10" 3 . Thus, our determination of 
2.0 ± 0.6 x 10" 3 is in excellent agreement with the literature de- 
terminations, and we have the advantage that we have estimated 
our stellar M/L ratios robustly, and our systematic uncertainties 
are better understood. We choose not to explore the quantita- 
tive consistency between our results and the integrated cosmic 
SFH, in part because the cosmic SFH is still poorly constrained 
in the particularly important epoch 0.5 < z < 2, where much of 
the stellar content in the Universe appears to have formed (e.g., 
Madau et al. 1996; Blain et al. 1999; Haarsma et al. 2000). 

Bell et al. (2003a) discuss implications of this result for the 
baryonic MF of galaxies, the cosmic mean density of stellar 
and cold gas in the local Universe, and the mean 'cold gas 
fraction' in the local Universe, accounting again for the main 
sources of systematic uncertainty. Salucci & Persic (1999) 
present a complementary analysis of spiral galaxies using pri- 
marily dynamically-derived M/L ratios, finding a similar but 
slightly higher cosmic mean density of stellar and cold gas mass 
compared to Bell et al. (2003a). Either way, both works show 
the exciting potential offered by a detailed understanding of the 
mass-to-light ratios and stellar MFs in the local Universe. 

It is interesting to use the stellar mass density estimate in 
conjunction with the luminosity densities that we derive to es- 
timate the stellar M/L ratio of the local Universe. Taking the 
stellar mass density of the Universe to be 5.5 ± 0.1 ± 1.6 x 
lO 8 /zM Mpc -3 , we find (3.64,3.50,3.05,2.57,2.00,0.95) ± 
0.03M Q /L Q in ugrizK passbands, with systematic errors of 
~ 30%, owing to uncertainties in SFHs and dust. These val- 
ues agree well with 0.93M Q /L Q from Cole et al. (2001) using 
2MASS TT-band (although we account for systematic error in 
the stellar M/L ratios), and the range 2.6-5. 2Mq/Lq estimate 
of Glazebrook et al. (2003) in the r-band using the cosmic mean 
galaxy spectrum as a constraint. Note that we scale both cases 
down from a Salpeter IMF to our maximally-massive 'diet' 
Salpeter IMF. Glazebrook et al. (2003) find a slightly larger 
range of possible M/L ratios than we do. Aperture effects may 
play a role, since they fit SFHs to the 3"fiber cosmic mean spec- 
trum, which will be biased towards the inner (redder) parts of 
galaxies (e.g., Bell & de Jong 2000). Another systematic dif- 
ference is that Glazebrook et al. construct the M/L ratio from 
the cosmic mean spectrum, which increases the uncertainties 
in how much light comes from old stellar populations relative 
to the younger stars with lower M/L ratio. We do not suffer 
as much from this uncertainty because we see much of the old 
stellar population light directly from red early-type galaxies. 12 

7. CONCLUSIONS 

A detailed understanding of the optical and NIR LFs and 
stellar MF is of fundamental importance to our understand- 
ing galaxy formation and evolution. We use a large sample of 
galaxies from 2MASS and the SDSS EDR to estimate the op- 
tical and NIR LFs and stellar MF in the local Universe, assum- 
ing a universally-applicable stellar IMF. We correct 2MASS ex- 
tended source catalog Kron magnitudes to total by comparing 

12 Although note that Rudnick et al. (2003) find that the sum of individual 
color-derived stellar masses for intensely star-forming galaxies can overpredict 
the total stellar mass by nearly a factor of two, whereas estimating the stellar 
mass from the summed light of these star-forming galaxies yields a much more 
accurate estimate, as the composite SF history of many bursty galaxies is more 
nearly smooth than the SF histories of each individual galaxy. This uncertainty 
will affect us only minimally, bearing in mind that the total stellar mass is domi- 
nated by early-type galaxies at the present epoch (>50%), and that the fraction 
of galaxies which are actively star-bursting is very low at the present epoch 
(e.g., Wolf etal. 2003). 
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2MASS magnitudes with deeper data from a variety of sources. 
Using the complete sky coverage of 2MASS to maximal effect, 
we find that the SDSS EDR region is <~ 8% overdense, and we 
correct for this overdensity when deriving LFs and stellar MFs. 

We estimate ^-corrections, evolution corrections, and present- 
day stellar M/L ratios by comparison of galaxy ugrizK magni- 
tudes with galaxy evolution models at each galaxy's redshift. 
The corrections and M/L ratios we derive are in excellent agree- 
ment with previous work, having the advantage that they incor- 
porate the multi-passband information to maximal effect. We 
estimate ^25% systematic errors in stellar M/L ratios from the 
effects of dust and bursts of SF, and we incorporate other ran- 
dom sources of error in our analysis. Assuming a universally- 
applicable stellar IMF, the limiting uncertainty in stellar M/L 
ratios is the overall uncertainty in the number of faint, low- 
mass stars. We conservatively adopt an IMF that has as many 
low-mass stars as possible without violating dynamical con- 
straints in nearby galaxies. Of course, IMFs poorer in low-mass 
stars are possible, and can be well-approximated by subtracting 
a logarithmic constant from all stellar M/L ratios and stellar 
masses presented in this paper. 

We construct optical and NIR LFs for galaxies in the local 
Universe using the V /V mix formalism, which has the key ad- 
vantage that it does not assume that galaxies in all environments 
have the same LF. The optical and NIR LFs that we estimate 
for this sample of galaxies agree to within the uncertainties 
with most recent literature optical and NIR LFs. We argue that 
2MASS misses faint, LSB galaxies, leading to underestimates 
of up to 25% in the TT-band luminosity density. The optical and 
NIR luminosity densities in the local Universe look to be well- 
constrained to within <~ 20%, and matches qualitatively a 'cos- 
mic SFH' that peaks at early times and continues to the present 
day at a reduced rate, with metallicities of roughly solar. 

We estimate the stellar MF using both TT-band and g-band 
limited galaxy samples. The MFs we derive using both sam- 
ples are consistent, given the bias against faint LSB galaxies in 
the 2MASS catalog. The g-band derived stellar MF goes down 
to lower stellar masses than the TT-band MF. Both MFs agree 
with the Cole et al. (2001) stellar MF, if the same IMFs are as- 
sumed. The faint end slope of the stellar MF is steeper than 
a = -1 . 1, reflecting the relatively low stellar M/L ratios of low 
mass galaxies. 13 

Under the assumption of a universally-applicable stellar IMF, 
the stellar mass density in the local Universe is Q^li = 2.0 ± 
0.6 x 10~ 3 , accounting for all sources of systematic and ran- 
dom error, except for IMF uncertainty. To change to a dif- 
ferent IMF, we would reduce this estimate by a constant fac- 
tor, as given in the text. Our stellar mass density estimate 
is consistent with earlier estimates, with the advantage that 
the systematic trends and errors in stellar M/L ratios are ac- 
counted for in this work. We find 'cosmic' stellar M/L ratios of 
(3.64, 3.50, 3.05, 2.57, 2.00, 0.95), in solar units in ugrizK, with 
~ 30% systematic errors, owing to the possible effects of bursts 
of SF and dust. 

Finally, we examine type-dependence in the optical and NIR 
LFs and the stellar MF. In agreement with previous work, we 
find that the characteristic luminosity or mass of early-type 
galaxies is larger than for later types, and the faint end slope 

13 Catalogs of stellar, estimated Hi and H2 gas masses (following 
Bell et al. 2003a) for 13 < r < 17.5 SDSS EDR galaxies, along 
with electronic versions of the tables in this paper, are given at: 

http : //www . mpia-hd . mpg . de/ homes /bell /data .html 



is steeper for later types than for earlier types. These results are 
robust to systematic differences in galaxy typing, although the 
overall numbers of early and late-type galaxies are somewhat 
dependent on the exact typing algorithm. Accounting for typ- 
ing uncertainties, we estimate that at least half, and perhaps as 
much as 3/4, of the stellar mass in the Universe is in early-type 
galaxies. 
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FIG. A19. — Comparison of the real g-band (open circles) and K-band (solid 
circles) LF, in black, with the predicted LFs transformed from the stellar MF 
using the average stellar M/L ratio at a given stellar mass (in grey). 
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APPENDIX 

MASS-TO-LIGHT RATIOS AND GALAXY COLORS 

As an aid to workers in the field, we present the K and g- 
band distributions of stellar M/L ratio, and the color-M/L ratio 
relations in the different SDSS and 2MASS passbands. 

The Distribution of Stellar M/L ratios 

In Fig. A18, we show the number density of galaxies as 
a function of their estimated stellar M/L ratio, assuming our 
maximal 'diet' Salpeter IMF. We denote different mass bins (in 
Mq) with different line styles: dot-dashed 9 < log 10 M/z 2 < 9.5, 
solid 9.5 < \og lQ Mh 2 < 10, dotted 10 < \og m Mh 2 < 10.5, and 
dashed 10.5 < log w Mh 2 < 11. In K-bznd, the 9 < \og m Mh 2 < 
9.5 bin is missing because of poor number statistics. In g-band, 
it is clear that the average stellar M/L ratio increases with in- 
creasing galaxy stellar mass, which indicates that more of the 
stars were formed at an earlier time (e.g., Kauffmann et al. 



Galaxy Luminosities and Stellar Masses 



21 




-0.4 -0.2 0.0 0.2 -0.5 0.0 0.5 1.0 

log 10 (M/L K ) log 10 (M/L g ) 

FIG. Al 8. — Distributions of stellar M/L ratio estimated from galaxy colors in ^-band (left) and g-band (right). We show four different galaxy stellar mass bins in 
units of solar mass (M Q ): 9 < log 10 M/i 2 < 9.5 (dot-dashed), 9.5 < log 10 M/i 2 < 10 (solid), 10 < log w Mh 2 < 10.5 (dotted), and 10.5 < log 10 M/i 2 < 11 (dashed). 
The A'-band 9 < log lQ Mh 2 < 9.5 bin is missing owing to poor number statistics. 



2003b). Moreover, the scatter becomes somewhat narrower in 
stellar M/L ratio at high stellar mass indicating less diversity in 
SFH. Massive galaxies tend to be rather old, regardless of mor- 
phological type, whereas less massive galaxies can have a wide 
range of ages (see also, e.g., Fig. 2 of Kauffmann et al. 2003b). 
Hints of this trend in /f-band are visible, but much weaker, 
showing the well-documented lack of sensitivity of M/Lk ra- 
tio to SFH (Bell & de Jong 2001). Sources of error include 
uncertainties from magnitude errors, systematic uncertainties 
in stellar M/L ratio from dust and bursts of SF (~ 25% in terms 
of stellar M/L ratio), and Poisson errors. The systematic un- 
certainties dominate, but are difficult to meaningfully estimate; 
thus, error bars are not given in this particular case. We tabulate 
the g and TT-band distributions in Table A6. 

It is worth briefly commenting on why these distributions 
are useful. In Fig. A19, we show the observed g-band (open 
circles) and /T-band (filled circles) LFs in black. Overplot- 
ted in grey are the predictions from the ^-band-derived stellar 
MF, using the average stellar M/L ratio as a function of stellar 
mass (as is often used by galaxy modelers to transform a stellar 
mass distribution into a luminosity function 14 ). To transform 
the stellar mass into luminosities, we adopt the bi-weight fit of 
stellar M/L ratio as a function of stellar mass: log 10 (M/Lp = 

-2.61 +0.298 log lo (M*/z 2 /M ), and log^M/L^) = -0.42 
+0.033 log w (M*h 2 /M e ). The estimated /T-band LF is in excel- 
lent agreement with the observed LF around the knee of the LF, 
and is in reasonable agreement at all luminosities with the pre- 
dicted /T-band LF, once 2MASS's selection bias against LSB 
galaxies is corrected for (the thick grey dashed line in Fig. 8). 
This indicates that the variation in stellar M/L ratio in TT-band at 
a given stellar mass is sufficiently small so that the predicted LF 
is close to the real LF. In contrast, using the average M/L g ratio 
at a given stellar mass is clearly insufficient to reproduce the g- 
band LF; especially at the faint end where one sees in Fig. A18 

14 Often, a single 'typical' stellar M/L ratio is used, which is even worse than 
the case we explore. 



that the scatter in M/L ? ratio is particularly large. This shows 
the importance of accounting for the spread in stellar M/L ratio 
at a given mass when transforming a stellar mass function into 
a luminosity function, especially in the optical. 15 

The Stellar M/L Ratio-Color Correlation 

Bell & de Jong (2001) presented relationships between John- 
son/Cousins optical-NIR colors and stellar M/L ratios in the 
optical and NIR using the galaxy models of Bell & Bower 
(2000). In this paper, we construct stellar M/L ratio estimates 
using galaxy evolution model fits to SDSS ugriz and 2MASS 
TT-band fluxes, finding tight correlations between optical color 
and stellar M/L ratio (e.g. Fig. 6). Therefore, in this appendix 
we compare our results with Bell & de Jong (2001), we exam- 
ine the color-M/L ratio correlations in detail where necessary, 
and we present fits to the color-M/L ratio correlations in the 
SDSS/2MASS passbands. 

To facilitate inter-comparison between our results, we choose 
to predict Johnson/Cousins M/L ratio values for the best-fit 
SEDs to the SDSS/2MASS galaxies in the g-band selected 
sample. We estimate B-R colors using the transformation 
in Fukugita et al. (1996): B-R = 1.506(g-r) + 0.370, with a 
~ 0.05 mag systematic error. 

We show two representative color-M/L ratio correlations in 
Fig. A20, where we show the M/L ratio in B-band (left-hand 
panel) and TT-band (right-hand panel) as a function of B-R 
color. Additionally, we give the least-squares 'robust' bi-weight 
fit to the estimates (solid line), and the relationship from the 
mass-dependent formation epoch with bursts model from Bell 
& de Jong (2001, dashed line). These two panels are easily 
compared with panel (d) of Fig. 1 in Bell & de Jong (2001); 
indeed, this is why we chose to estimate B-R colors and John- 
son/Cousins M/L ratios. The agreement between the ieast- 

I5 This also applies to transforming LFs from one passband into another using 
a luminosity-dependent typical color. If the spread in color at a given luminosity 
is rather small, one can get away with this approximation. If the spread in color 
is larger, the estimated LF will be biased. 




Table A7 

Stellar M/L ratio as a function of color 



Color 


Clg 


h 


Clr 


br 


a; 


b, 




b z 


aj 


bj 


a H 


b H 


a K 


b K 


u-g 


-0.221 


0.485 


-0.099 


0.345 


-0.053 


0.268 


-0.105 


0.226 


-0.128 


0.169 


-0.209 


0.133 


-0.260 


0.123 


u — r 


-0.390 


0.417 


-0.223 


0.299 


-0.151 


0.233 


-0.178 


0.192 


-0.172 


0.138 


-0.237 


0.104 


-0.273 


0.091 


u-i 


-0.375 


0.359 


-0.212 


0.257 


-0.144 


0.201 


-0.171 


0.165 


-0.169 


0.119 


-0.233 


0.090 


-0.267 


0.077 


u-z 


-0.400 


0.332 


-0.232 


0.239 


-0.161 


0.187 


-0.179 


0.151 


-0.163 


0.105 


-0.205 


0.071 


-0.232 


0.056 


g-r 


-0.499 


1.519 


-0.306 


1.097 


-0.222 


0.864 


-0.223 


0.689 


-0.172 


0.444 


-0.189 


0.266 


-0.209 


0.197 


g-i 


-0.379 


0.914 


-0.220 


0.661 


-0.152 


0.518 


-0.175 


0.421 


-0.153 


0.283 


-0.186 


0.179 


-0.211 


0.137 


g-z 


-0.367 


0.698 


-0.215 


0.508 


-0.153 


0.402 


-0.171 


0.322 


-0.097 


0.175 


-0.117 


0.083 


-0.138 


0.047 


r-i 


-0.106 


1.982 


-0.022 


1.431 


0.006 


1.114 


-0.052 


0.923 


-0.079 


0.650 


-0.148 


0.437 


-0.186 


0.349 


r-z 


-0.124 


1.067 


-0.041 


0.780 


-0.018 


0.623 


-0.041 


0.463 


-0.011 


0.224 


-0.059 


0.076 


-0.092 


0.019 


Color 


a B 


b B 


av 


bv 


Or 


b R 


at 


bi 


aj 


bj 


a H 


b H 


a K 


b K 


B-V 


-0.942 


1.737 


-0.628 


1.305 


-0.520 


1.094 


-0.399 


0.824 


-0.261 


0.433 


-0.209 


0.210 


-0.206 


0.135 


B-R 


-0.976 


1.111 


-0.633 


0.816 


-0.523 


0.683 


-0.405 


0.518 


-0.289 


0.297 


-0.262 


0.180 


-0.264 


0.138 



Note. — Stellar M/L ratios are given by log 10 (M/L) = a\ + {b\ X Color) where the M/L ratio is in solar units. If all galaxies are sub-maximal then the above zero 
points (a\) should be modified by subtracting an IMF dependent constant as follows: 0.15 dex for a Kennicutt or Kroupa IMF, and 0.4 dex for a Bottema IMF. Scatter 
in the above correlations is ~ 0.1 dex for all optical M/L ratios, and 0.1-0.2 dex for NIR M/L ratios (larger for galaxies with blue optical colors). SDSS filters are in 
AB magnitudes; Johnson BVR and JHK are in Vega magnitudes. 



squares fit to many passbands' methodology of this paper is 
in excellent agreement with the galaxy modeling of Bell & de 
Jong (2001) for the B-band M/L ratio-color relation. This re- 
sult is insensitive to detailed passband choice since all the op- 
tical M/L ratio-color correlations that we derive are consistent 
with Bell & de Jong (2001). 

In the right-hand panel of Fig. A20, we show the run of K- 
band stellar M/L ratio estimates from the 'least-squares fit to 
many passbands' methodology we adopt in this paper (points) 
against the galaxy model-based estimate of Bell & de Jong 
(2001, dashed line). We note the somewhat poorer agree- 
ment between these two different methodologies. There is a 
zero point offset, owing to our use of PEGASE rather than the 
Bruzual & Chariot (in preparation) models. Furthermore, there 
is considerably more scatter at the blue end of the correlation 
than the models of Bell & de Jong (2001) predict, and a some- 
what shallower correlation than is expected on the basis of their 
galaxy modeling. The data points fill in the range of possible 
colors and M/L ratios of stellar populations with a wide variety 
of ages and metallicities Bell & de Jong (2001, Fig. 2, panel 



c), indicating that this spread is primarily caused by a spread in 
metallicity . In particular, it is clear that optically-blue galax- 
ies have a wide range of estimated metallicities. Recall that 
our methodology in this paper is to estimate ages and metal- 
licities using the optical-NIR colors following Bell & de Jong 
(2000). In contrast, the galaxy evolution models of Bell & de 
Jong (2001) do not have a large metallicity spread, and there- 
fore, do not reproduce this feature. This large metallicity spread 
flattens the optical color-M/L/j ratio correlation, adding ~ 0.2 
dex scatter at the blue end. It is important to note that this 0.2 
dex scatter is no more than a factor of two in excess of the scat- 
ter in optical M/L ratios as a function of color. 

In conclusion, we find that the optical M/L ratios as a func- 
tion of color are in good agreement with Bell & de Jong (2001). 
However, for NIR M/L ratios we find that real galaxies sug- 
gest a larger metallicity scatter than accounted for by Bell & de 
Jong (2001), leading to a shallower color- M/L ratio slope and a 
larger spread at the blue end. In Table A7 we present the cor- 
relations between SDSS ugriz colors and SDSS/2MASS M/L 
ratios, and between BVR colors and Johnson/Cousins M/L ra- 
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tios, to allow intercomparison with Bell & de Jong (2001), and 
to allow the estimation of systematic differences between their 
work and ours in the NIR. Typical M/L ratio uncertainties are 
^0.1 dex in the optical, and 0.1 (0.2) dex in the NIR at the red 
(blue) end. We do not present M-band M/L estimates for M-band 
(because of its strong sensitivity to recent SF) or correlations 
involving only NIR colors (because of their strong metallicity 
sensitivity). 

Table A6 
Stellar M/L ratio distributions 



M/L 


log lo (0//! J Mpc- 


" J log 10 M- 


nog l0 [M/LT L ) 


(1) 


(2) 


(3) 


(4) 


(5) 


K-band 


-0.40 






-3.65 




-0.36 




-3.08 






-0.32 


-2.16 




-3.21 




-0.28 






-3.39 


-4.35 


-0.24 


-1.64 


-2.16 


-2.69 


-3.53 


-0.20 


-1.97 


-1.69 


-2.11 


-2.90 


-0.16 


-1.74 


-1.45 


-1.67 


-2.24 


-0.12 


-1.83 


-1.46 


-1.31 


-1.72 


-0.08 


-1.43 


-1.34 


-1.26 


-1.50 


-0.04 


-1.86 


-1.41 


-1.29 


-1.50 


0.00 


-2.16 


-1.72 


-1.66 


-1.68 


0.04 


-1.66 


-1.79 


-2.21 


-2.67 


0.08 


-1.49 


-1.83 


-2.82 


-3.36 


0.12 


-1.50 


-2.17 


-3.45 


-3.58 


0.16 




-2.12 






g-band 


-0.46 


-4.07 








-0.37 


-3.29 








-0.29 


-3.34 


-4.05 






-0.20 


-2.73 


-3.50 


-4.57 




-0.12 


-2.18 


-2.97 


-4.25 




-0.03 


-1.89 


-2.48 


-3.53 




0.05 


-1.66 


-2.14 


-3.06 


-4.96 


0.14 


-1.49 


-1.92 


-2.64 


-4.41 


0.22 


-1.57 


-1.81 


-2.31 


-3.88 


0.31 


-1.53 


-1.77 


-2.12 


-3.02 


0.39 


-1.60 


-2.04 


-2.15 


-2.78 


0.48 


-1.62 


-1.71 


-1.94 


-2.45 


0.56 


-1.80 


-1.60 


-1.65 


-2.06 


0.65 


-2.25 


-1.83 


-1.56 


-1.77 


0.73 


-2.20 


-2.48 


-1.97 


-1.87 


0.82 




-3.09 


-2.58 


-2.44 


0.90 




-2.92 


-3.23 


-3.43 


0.99 


-1.81 


-2.13 


-2.58 


-3.32 



Note. — The distribution of stellar M/L ratio (1) in solar units is given 
for different bins of mass in solar units: 
(2) 9 < \og l0 Mh 2 < 9.5; (3) 9.5 < log 10 M/i 2 < 10; 
(4) 10 < log 10 Mh 2 < 10.5; (5) 10.5 < log 10 M/i 2 < 11 



